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Abstract 

We consider many-fermion systems with singular Fermi surfaces, which 
contain Van Hove points where the gradient of the band function k i— > e(k) 
vanishes. In a previous paper, we have treated the case of spatial dimen- 
sion d > 3. In this paper, we focus on the more singular case d = 2 and 
establish properties of the fermionic self-energy to all orders in perturba- 
tion theory. We show that there is an asymmetry between the spatial and 
frequency derivatives of the self-energy. The derivative with respect to the 
Matsubara frequency diverges at the Van Hove points, but, surprisingly, the 
self-energy is C 1 in the spatial momentum to all orders in perturbation the- 
ory, provided the Fermi surface is curved away from the Van Hove points. 
In a prototypical example, the second spatial derivative behaves similarly to 
the first frequency derivative. We discuss the physical significance of these 
findings. 



*feldman@math.ubc.ca; supported by NSERC of Canada 

tsalmhofer@itp.uni-leipzig.de; supported by DFG grant Sa-1362, an ESI senior research fel- 
lowship, and NSERC of Canada 



1 



Contents 



1 Introduction 3 

2 Main Results 5 

2. 1 Hypotheses on the dispersion relation 6 

2.2 Main theorem 7 

2.3 Heuristic explanation of the asymmetry 8 

3 Fermi Surface 11 

3.1 Normal form for e(k) near a singular point 11 

3.2 Length of overlap estimates 13 

3.2.1 Length of overlap - special case 14 

3.2.2 Length of overlap - general case 16 

4 Regularity 23 

4.1 The gradient of the self-energy 23 

4.1.1 The second order contribution 23 

4.1.2 The general diagram 29 

4.2 The frequency derivative of the self-energy 35 

5 Singularities 38 

5.1 Preparations 39 

5.2 g -derivative 41 

5.3 First spatial derivatives 45 

5.4 The second spatial derivatives 46 

5.5 One-loop integrals for the xy case 58 

5.5.1 The particle-hole bubble 58 

5.5.2 The particle-particle bubble 59 

6 Interpretation 60 

6.1 Asymmetry and Fermi velocity suppression 62 

6.2 Inversion problem 62 

A Interval Lemma 65 

B Signs etc. 66 



2 



1 Introduction 



In this paper, we continue our analysis of (all-order) perturbative properties of 
the self-energy and the correlation functions in fermionic systems with a fixed 
non-nested singular Fermi surface. That is, the Fermi surface contains Van Hove 
points, where the gradient of the dispersion function vanishes, but satisfies a no- 
nesting condition away from these points. In a previous paper [1], we treated the 
case of spatial dimensions d > 3. Here we focus on the two-dimensional case, 
where the effects of the Van Hove points are strongest. 

We have given a general introduction to the problem and some of the main 
questions in [1]. As discussed in [1], the no-nesting hypothesis is natural from 
a theoretical point of view, because it separates effects coming from the saddle 
points and nesting effects. Moreover, generically, nesting and Van Hove effects 
do not occur at the same Fermi level. In the following, we discuss those aspects 
of the problem that are specific to two dimensions. As already discussed in [1], 
the effects caused by saddle points of the dispersion function lying on the Fermi 
surface are believed to be strongest in two dimensions (we follow the usual jargon 
of calling the level set the Fermi "surface" even though it is a curve in d = 2). 
Certainly, the Van Hove singularities in the density of states of the noninteracting 
system are strongest in d = 2. As concerns many-body properties, we have shown 
in [1] that for d > 3, the overlapping loop estimates of [10] carry over essentially 
unchanged, which implies differentiability of the self-energy and hence a quasi- 
particle weight (Z-factor) close to 1 to all orders in renormalized perturbation 
theory. 

In this paper, we show that for d = 2, there are more drastic changes. Namely, 
there is an asymmetry between the derivatives of the self-energy S(g , q) with 
respect to the frequency variable q and the spatial momentum q. We prove that 
the spatial gradient VS is a bounded function to all orders in perturbation theory 
if the Fermi surface satisfies a no-nesting condition. By explicit calculation, we 
show that for a standard saddle point singularity, even the second-order contribu- 
tion do^iqo, q s ) diverges as (log |go|) 2 at any Van Hove point q s (if that point is 
on the Fermi surface). 

This asymmetric behaviour is unlike the behaviour in all other cases that are 
under mathematical control: in one dimension, both <9 £2 and <9iE 2 diverge like 
log \ q \ at the Fermi point. This is the first indication for vanishing of the Z- 
factor and the occurrence of anomalous decay exponents in this model. The point 
is, however, that once a suitable Z-factor is extracted, Zd{L 2 remains of order 
1 in one dimension, while for two-dimensional singular Fermi surfaces, the p- 



3 



dependent function Z(p)VE(p) vanishes at the Van Hove points. In higher di- 
mensions d > 2, and with a regular Fermi surface fulfilling a no-nesting condition 
very similar to that required here, £ is continuously differentiable both in q and 
in q. Thus it is really the Van Hove points on the Fermi surface that are responsi- 
ble for the asymmetry. In the last section of this paper, we point out some possible 
(but as yet unproven) consequences of this behaviour. 

Our analysis is partly motivated by the two-dimensional Hubbard model, a 
lattice fermion model with a local interaction and a dispersion relation k i— > e(k) 
which, in suitable energy units, reads 

e(k) = — cos ki — cos k 2 + 9(1 + cos k x cos k 2 ) — /i. (1) 

The parameter fi is the chemical potential, used to adjust the particle density, and 
9 is a ratio of hopping parameters. As we shall explain now, the most interesting 
parameter range is /x ps and 

< 9 < 1. 

The zeroes of the gradient of e are at (0, 0), (tt, tt) and at (tt, 0), (0, tt). The first 
two are extrema, and the last two are the saddle points relevant for van Hove 
singularities (VHS). For = 0, both saddle points are on the Fermi surface. For 
9 = 1 the Fermi surface degenerates to the pair of lines {ki = 0} U {k 2 = 0}, 
so we assume that 9 < 1. For 9 = and /i = 0, the Fermi surface becomes the 
so-called Umklapp surface U = {k : ki ± k 2 = ±7r}, which is nested since it 
has flat sides. For < 9 < 1 the Fermi surface at fx = has nonzero curvature 
away from the Van Hove points (tt, 0) and (0, tt). Viewed from the point (it, tt), it 
encloses a strictly convex region (as a subset of M 2 ). There is ample evidence that 
in the Hubbard model, it is the parameter range 9 > and electron density near to 
the van Hove density (p, ps 0) that is relevant for high-T c superconductivity (see, 
e.g. [2, 3, 4, 5, 6]). In this parameter region, an important kinematic property is 
that the two saddle points at (tt, 0) and (0, tt) are connected by the vector Q = 
(tt, tt), which has the property that 2Q = mod IttI?. This modifies the leading 
order flow of the four-point function strongly (Umklapp scattering, [4, 5, 6]). The 
bounds we discuss here hold both in presence and absence of Umklapp scattering. 

The interaction of the fermions is given by \v, where A is the coupling constant 
and v is the Fourier transform of the two-body potential defining a density-density 
interaction. For the special case of the Hubbard model, two fermions interact only 
if they are at the same lattice point, so that v(k) = 1. Despite the simplicity of the 
Hamiltonian, little is known rigorously about the low-temperature phase diagram 
of the Hubbard model, even for small |A|. In this paper, we do perturbation theory 
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to all orders, i.e. we treat A as a formal expansion parameter. For a discussion 
of the relation of perturbation theory to all orders to renormalization group flows 
obtained from truncations of the RG hierarchies, see the Introduction of [1]. 

Although our analysis is motivated by the Hubbard model, it applies to a much 
more general class of models. In this paper, we shall need only that the band 
function e has enough derivatives, as stated below, and a similar condition on 
the interaction. In fact, the interaction is allowed to be more general than just a 
density-density interaction: it may depend on frequencies, as well as the spin of 
the particles. See [10, 11, 12, 13] for details. As far as the singular points of e are 
concerned, we require that they are nondegenerate. The precise assumptions on e 
will be stated in detail below. 

We add a few remarks to put these assumptions into perspective. No matter if 
we start with a lattice model or a periodic Schrodinger operator describing Bloch 
electrons in a crystal potential, the band function given by the Hamiltonian for 
the one-body problem is, under very mild conditions, a smooth, even analytic 
function. In such a class of functions, the occurrence of degenerate critical points 
is nongeneric, i.e. measure zero. In other words, if 

e(k s + Rk) = -e x k\ + e 2 k\ + ... 

around a Van Hove point k s (here R is a rotation that diagonalizes the Hessian 
at k s ), getting even one of the two prefactors £; to vanish in a Taylor expansion 
requires a fine-tuning of the hopping parameters, in addition to the condition that 
the VH points are on S. Thus, in a one-body theory, an extended VHS, where the 
critical point becomes degenerate because, say, e\ vanishes, is nongeneric. On the 
other hand, experiments suggest [3] that e\ is very small in some materials, which 
seem to be modeled well by Hubbard-type band functions. On the theoretical side, 
in a renormalized expansion with counterterms, it is not the dispersion relation of 
the noninteracting system, but that of the interacting system, which appears in all 
fermionic co variances. It is thus an important theoretical question to decide what 
effects the interaction has on the dispersion relation and in particular whether an 
extended VHS can be caused by the interaction. We shall discuss this question 
further in Section 6. 

2 Main Results 

In this section we state our hypotheses on the dispersion function and the Fermi 
surface, and then state our main result. 
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2.1 Hypotheses on the dispersion relation 

We make the following hypotheses on the dispersion relation e and its Fermi sur- 
face T = { k | e(k) = } in d = 2. 

HI { k | |e(k)| < 1 } is compact. 
H2 e(k)isC r withr > 7. 

H3 e(k) = and Ve(k) = simultaneously only for finitely many k's, called 
Van Hove points or singular points. 

H4 If k is a singular point then [ f k dk e(k)] 1<j - <d is nonsingular and has one 
positive eigenvalue and one negative eigenvalue. 

H5 There is at worst polynomial flatness. This means the following. Let k E T . 
Suppose that k 2 — k 2 — /(&a — k\) is a C r ~ 2 curve contained in T in a 
neighbourhood of k. (If k is a singular point, there can be two such curves.) 
Then some derivative of f(x) at x — of order at least two and at most r — 2 
does not vanish. Similarly if the roles of the first and second coordinates are 
exchanged. 

H6 There is at worst polynomial nesting. This means the following. Let k e T 
and p G T with p. Suppose that k 2 — k 2 = f(ki — k\) is a C r ~ 2 curve 
contained in T in a neighbourhood of k and k 2 — p 2 = g(k\—pi) is a C r ~ 2 
curve contained in T in a neighbourhood of p. Then some derivative of 
f(x) — g(x) at x = of order at most r — 2 does not vanish. Similarly if the 
roles of the first and second coordinates are exchanged. If e(k) is not even, 
we further assume a similar nonvanishing when / gives a curve in T in a 
neighbourhood of any k e T and g gives a curve in — T in a neighbourhood 
of any p e — J 7 . 

We denote by n the largest nonflatness or nonnesting order plus one, and assume 
that 

r > 2n + 1 

The Fermi surface for the Hubbard model with < 9 < 1 and fx — 0, when 
viewed from (7r, 7r), encloses a convex region. See the figure below. It has nonzero 
curvature except at the singular points. If one writes the equation of (one branch 
of) the Fermi surface near the singular point (0, it) in the form k 2 — n = f(k{), 
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then /( 3 )(0) 7^ 0. So this Fermi surface satisfies the nonflatness and no-nesting 
conditions with n = 4. 



2.2 Main theorem 

In the following, we state our main results about the fermionic self-energy. A 
discussion will be given at the end of the paper, in Section 6. 

Theorem 2.1 Let B = R 2 /27rZ 2 and e G C 7 (B, R). Assume that the Fermi 
surface S = {k G B : e(k) = 0} contains points where Ve(k) = 0, and that 
the Hessian of 'e at these points is nonsingular. Moreover, assume that away from 
these points, the Fermi surface can have at most finite-order tangencies with its 
(possibly reflected) translates and is at most polynomially flat. [These hypotheses 
have been spelled out in detail in H1-H6 above.] As well, the interaction v is 
assumed to be short-range, so that its Fourier transform v is C 2 . 

Then there is a counterterm function K G C 1 (i3, R), given as a formal power 
series K = J2 r >i K r \ r in the coupling constant \, such that the renormalized 
expansion for all Green functions, at temperature zero, is finite to all orders in A. 

1. The self-energy is given as a formal power series E = Xl r >i ^r-A r , where 
for all r G N and all u G R, the function k i— > lu r {u),\z) G C X {B, C). 
Specifically, we have 



with the constants depending on r. Moreover, the function uj \— > S r (u;, k) is 
C 1 in uj for all k G B\V, where V denotes the integer lattice generated by 
all van Hove points, and the bar means the closure in B. 





< const 



< const 
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2. For e given by the normal form e(k) = kik 2 , which has a van Hove point at 
k = 0, the second order contribution S 2 to the self-energy obeys 

Zrod w £ 2 (c;,0) = -a x (logM) 2 + 0(|logM|) 
Red fel d fe2 £ 2 (u;,0) = a 2 (log \u\f + 0(| log |w||) 

<9f £2(^5 0) anJ (9| 2 S2(w, 0) grow at most linearly in log The explicit 
values of a\ > and a 2 > are given in Lemmas 5.1 and 5.2, below. 

Theorem 2. 1 is a statement about the zero temperature limit of E. That is, £ r (uj, k) 
and its derivatives are computed at a positive temperature T = (3" 1 , where they 
are C 2 in uj and q, and then the limit (3 — > oo is taken. (Because only one- 
particle-irreducible graphs contribute to £ r , it is indeed a regular function of uj 
for all e R at any inverse temperature f3 < oo.) 

The bounds to all orders stated in item 1 of Theorem 2.1 generalize to low 
positive temperatures in an obvious way: the length-of-overlap estimates and the 
singularity analysis done below only use the spatial geometry of the Fermi surface 
for e, which is unaffected by the temperature. The other changes are merely to 
replace some derivatives with respect to frequency by finite differences, which 
only leads to trivial changes. 

Our explicit computation of the asymptotics in the model case of item 2 of 
Theorem 2.1 uses that several contributions to these derivatives vanish in the limit 
P — ■> oo, and that certain cancellations occur in the remaining terms. For this 
reason, the result stated in item 2 is a result at zero temperature. (In particular, the 
coefficients in the 0(log |u;|) terms are just numbers.) However, we do not expect 
any significant change in the asymptotics at low temperature and small uj to occur. 
That is, we expect the low-temperature asymptotics to contain only terms whose 
supremum over > n/f3 is at most of order (log/3) 2 , and the square of the 
logarithm to be present. 

2.3 Heuristic explanation of the asymmetry 

We refer to the different behaviour of VS (which is bounded to all orders) and 
c^X (which is log 2 -divergent in second order) as the asymmetry in the derivatives 
of S. In the case of a regular Fermi surface, no-nesting implies that S r is in C 1+s 
with a Holder exponent 5 that depends on the no-nesting assumption. A similar 
bound was shown in [1] for Fermi surfaces with singularities in d > 3 dimen- 
sions. In [1], we formulated a slight generalization of the no-nesting hypothe- 
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sis of [10], and again proved a volume improvement estimate, which implies the 
above-mentioned Holder continuity of the first derivatives. 

In the more special case of a regular Fermi surface with strictly positive cur- 
vature, we have given, in [11, 12], bounds on certain second derivatives of the 
self-energy with respect to momentum. We briefly review that discussion for the 
second-order contribution, to motivate why there is a difference between the spa- 
tial and the frequency derivatives. 

For simplicity, we assume a local interaction, and consider the infrared part of 
the two-loop contribution 

J(9o,q) = (C(uj 1 ,e(p 1 ))C( y U2,e(p 2 ))C(u,e)) 

where u = qo + viu>i + v 2 uj 2 , = ±1 and 

e = e(q + wipi + v 2 p 2 ) 

The angular brackets denote integration of pi and p 2 over the two-dimensional 
Brillouin zone and Matsubara summation of ui and lu 2 , over the set j|(2Z + 1). 
By infrared part we mean that the fermion propagators are of the form 

where U is a suitable cutoff function that is supported in a small, fixed neighbour- 
hood of zero. 

The third denominator depends on the external momentum (q , q) and deriva- 
tives with respect to the external momentum increase the power of that denomi- 
nator, which may lead to bad behaviour as (5 — > 00. The main idea why some 
derivatives behave better than expected by simple counting of powers (see [11]) 
is that in dimension two and higher, there are, in principle, enough integrations 
to make a change of variables so that e, e(pi) and e(p 2 ) all become integration 
variables. This puts all dependence on the external variable q into the Jacobian 
J of this change of variables. If J were C k with uniform bounds, /(go? q) would 
be C k in the spatial momentum q, and (by integration by parts) also in q . How- 
ever, J always has singularities, and the leading contributions to the derivatives 
of I(q , q) come from the vicinity of these singularities. It was proven in [11, 12] 
that if the Fermi surface is regular and has strictly positive curvature, these singu- 
larities of the Jacobian are harmless, provided derivatives are taken tangential to 
the Fermi surface. 
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To explain the change of variables, we first show it for the case without van 
Hove singularities and then discuss the changes required when van Hove singu- 
larities are present. 

In a neighborhood of the Fermi surface, we introduce coordinates p and 9, 
so that p = P(p,9). The coordinates are chosen such that p = e(p) and that 
P(p, 9 + 7r) is the antipode of P(p, 9). (We are assuming that the Fermi surface is 
strictly convex - see [11].) Doing this for p x and p 2 , with corresponding Jacobian 
J = det P', we have 



I(Qo, q) 



d9 1 d9 2 J{pi,9 1 )J(p 2 ,9 2 ) C(u,e) 



1,2 



where [F]i )2 now denotes multiplying F by C(u\, pi)C(u 2 , p 2 ) and integrating 
over pi and p 2 and summing over the frequencies. To remove the q-dependence 
from C(u,e), one now wants to change variables from 9i or 9 2 to e. This works 
except near points where 

de de 
W 1 = d9~ 2 = ' 

These equations determine the singularities of the Jacobian. The detailed analysis 
of their solutions is in [11]. Essentially, if one requires that the momenta p x and 
p 2 are on the FS, i.e. pi = p 2 = 0, that q = P(0, q ) is on the FS and that 
the sum q + i>iP(0, 9i) + v 2 P(0, 9 2 ) is on the FS, then the only solutions are 
9\ G {6* q , 9^ + 7r} and 9 2 e {# q , # q + n}. (The general case of momenta near to 
the Fermi surface is then treated by a deformation argument which requires that 
Ve 7^ and that the curvature be nonzero.) A detailed analysis of the singularity 
in J, in which strict convexity enters again, then implies that the self-energy is 
regular. 

The conditions needed for the above argument fail at the Van Hove points. But, 
introducing a partition of unity on the Fermi surface, they still hold away from 
the singular points. So the only contributions that may fail to have derivatives 
come from q + v±pi + v 2 p 2 in a small neighbourhood of the singular points. 
When a derivative with respect to go is taken, the integrand contains a factor of 
— i(icj — e)~ 2 . When a derivative with respect to q is taken, the integrand contains 
a factor Ve(q + f ipi + f2p 2 )(i^ — e) -2 . Because we are in a small neighbourhood 
of the singular point, the numerator, Ve(q + viPi + v 2 p 2 ), in the latter expression 
is small, and vanishes at the singular point. This suggests that the first derivative 
with respect to q may be better behaved than the first derivative with respect to 
q , as is indeed the case - see item 1 of Theorem 2.1. A second derivative with 
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respect to q may act on the numerator, Ve(q + v 1 p 1 + V2P2), and eliminate its 
zero. This suggests that the second derivative with respect to q behaves like the 
first derivative with respect to q , as is indeed the case - see item 2 of Theorem 
2.1. 

The above heuristic discussion is only to provide a motivation as to why the 
asymmetry in the derivatives of S 2 occurs. The proof does not make use of the 
idea of the change of variables to e, but rather of length-of-overlap estimates, 
which partially replace the overlapping loop estimates, away from the singular 
points. This allows us to show the convergence of the first q-derivative under 
conditions H1-H6, which are significantly weaker than strict convexity, and it 
also allows us to treat the situation with Umklapp scattering, which had to be 
excluded in second order in [1 1], and which is the reason for the restriction on the 
density in [11]. 

3 Fermi Surface 

In this section we prove bounds on the size of the overlap of the Fermi surface 
with translates of a tubular neighbourhood of the Fermi surface. These bounds 
make precise the geometrical idea that for non-nested surfaces (here: curves), the 
non-flatness condition H5 strongly restricts such lengths of overlap. 

3.1 Normal form for e(k) near a singular point 

Lemma 3.1 Let d = 2 and assume H2-H5 with r > n + 1. Assume that k = 
is a singular point of e. Then there are 

o integers 2 < v 1 , u 2 < n , 

o a constant, nonsingular matrix A and 

o (7r-2-max{i>i,i/ 2 } f unC fi ons a(k), b(k) and c(k) that are bounded and bounded 
away from zero 

such that in a neighbourhood of the origin 

e{Ak) = a(k)(A; 1 - k^b{k)){k 2 - £fc(k)). (2) 
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Proof: Let Ai and A 2 be the eigenvalues of Qk , e(0)] 1<i j<2 and set A = 

a/|A 2 /Ai|. By the Morse lemma [14, Lemma 1.1 of Chapter 6], there is a C r ~ 2 
diffeomorphism x(k) with x(0) = such that 

e(k) = A 1 x 1 (k) 2 + A 2 a; 2 (k) 2 

= A 1 (x 1 (k) 2 -A 2 a : 2 (k) 2 ) 

= A 1 (x 1 (k)-Ax 2 (k))(x 1 (k) + Ax 2 (k)) 

= Ai(aiA;i + a 2 k 2 - x 1 (k) ) (hki + b 2 k 2 - x 2 (k) ) 

with xi(k) and x 2 (k) vanishing to at least order two at k = 0. Here aiki + a 2 k 2 
and biki + b 2 k 2 are the degree one parts of the Taylor expansions of Xi(k) — 
Xx 2 (k) and Xx(k) + Ax 2 (k) respectively. Since the Jacobian det Dx(0) of the 

diffeomorphism at the origin is nonzero and det 



A 



2A 7^ 0, we have 



det 



ai a 2 



det 



1 -A 
1 A 



L>x(0) + 



Setting 

we have 
Write 
with 



a\ a 2 
h b 2 



e(Ak) = X 1 (k 1 - xi(.4k))(A;2 - x 2 (Ak)) 
h - x t (Ak) = h - hfrQt) - k? 9l (k 2 ) 

kMk) = x^Ak) - ^(Ak)^ 
Wgifa) = 5i(Ak)| fci=0 

Since x\ (Ak) vanishes to order at least two at k = 0, 



Jo 



k=(tfci,A;2) 



dt 



is C r 3 and vanishes to order at least one at k = and, in particular, |/i(k) | < \ 
for all k in a neighbourhood of the origin. We choose v x to be the power of 
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the first nonvanishing term in the Taylor expansion of xi(Ak)\ ki=0 . Since this 
function must vanish to order at least two in k 2 , we have that v\ > 2. By the 
nonflatness condition, applied to the curve implicitly determined by k\ = x\ (Ak), 
v\ < r — 2. So gi(k 2 ) is C r ~ 2 ~ ui and is bounded and bounded away from zero in 
a neighbourhood of k 2 = 0. 
In a similar fashion, write 

k 2 - x 2 (Ak) =k 2 - k 2 f 2 (k) - k\ 2 g 2 (ki) 

with 

k 2 f 2 (k) = x 2 (Ak) - x 2 (Ak)\ k2=0 
k? 92(h) = x 2 {Ak)\ k2=0 
Then we have the desired decomposition (2) with 

a(k) = Ai(l-/i(k))(l-/ 2 (k)) 

6(k) = 9l{k £ j . , c(k)= 92{h) 



l-/i(k) v ' l-/ 2 (k) 



We remark that it is possible to impose weaker regularity hypotheses by ex- 
ploiting that k v 2 b{ki), resp. fc^c(k), is a C r ~ 3 function whose /c 2 , resp. k\, deriva- 
tives of order strictly less than v x , resp. u 2 , vanish at k 2 = 0, resp. k x = 0. 



3.2 Length of overlap estimates 

It follows from the normal form derived in Lemma 3.1 that under the hypotheses 
H2-H5 the curvature of the Fermi surface may vanish as one approaches the sin- 
gular points. Thus, even if the Fermi surface is curved away from these points, 
there is no uniform lower bound on the curvature. Curvature effects are very im- 
portant in the analysis of regularity estimates, and in a situation without uniform 
bounds these curvature effects improve power counting only at scales lower than 
a scale set by the rate at which the curvature vanishes. Thus it becomes natural to 
define, at a given scale, scale-dependent neighbourhoods of the singular points, 
outside of which curvature improvements hold. The estimates for the length of 
overlaps that we prove in this section allow us to make this idea precise. They 
hold under much more general conditions than a nonvanishing curvature, namely 
the nonnesting and nonflatness assumptions H5 and H6 suffice. We first discuss 
the special case corresponding to the normal form in the vicinity of a singular 
point, and then deal with the general case. 
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3.2.1 Length of overlap - special case 
Lemma 3.2 Let v\>2 and u 2 > 2 be integers and 

e(x,y) = {x- y ui b(x,y))(y - x U2 c(x,y)) 

with b and c bounded and bounded away from zero and with b,c G C U2+1 . Let 
u{x) obey 

u{x) = x V2 c(x,u(x)) 

for all x in a neighbourhood of 0. That is, y = u(x) lies on the Fermi curve 
e(x,y) = 0. There are constants C and D > such that for all e > and 
0<5< \(X,Y)\ < D 

Vo\{xeR \ \x\ <D, \e(X + x,Y + u(x))\ <e} < C(f ) 1/u2 

Proof: Write 

e(X + x, Y + u(x)) = F(x, X, Y)G{x, X, Y) (3) 

with 

F(x,X,Y) = X + x- (Y + u(x)Y 1 b(X + x,Y + u(x)) 
G(x,X,Y) = Y + u(x) - (X + x) V2 c(X + x,Y + u(x)) 
= Y- {(X + xY 2 - x U2 }c(X + x,Y + u(x)) 
—x V2 {c(X + x, Y + u(x)) — c(x, u(x))} 

Observe that, for all allowed x, X and Y, 

\F(x,X,Y)\<^5 \k F (x,X,Y)\>Z 

since x, X, Y and u(x) all have to be 0(D) small. For our analysis of G(x, X, Y) 
we consider two separate cases. 

Case 1: \Y\ > k\X\ with k a constant to be chosen shortly. Since c(X + x, Y + 
u(x)) — c(x, u(x)) vanishes to first order in (X, Y), for all x 

\x U2 {c(X + x,Y + u(x)) -c(x,u(x))}\ < y^[|X| + |F|] 
\l- x [x» 2 {c(X + x,Y + u(x))-c(x,u(x))}}\ < ^[\X\ + \Y\] 

Since (X + x) U2 — x u ' 2 vanishes to first order in X, for all x 

\{(X + xf 2 -x U2 }c(X + x,Y + u(x))\ < k\X\ 
\^[{(X + x) U2 -x U2 }c(X + x,Y + u(x))}\ < k\X\ 
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We choose k = max{2, 200/t}. Then 

\G(x,X,Y)\ > *L\Y\ \l- x G(x,X,Y)\ < ±\Y\ 

Thus, by (3) and the product rule, 

\&{X + x,Y + u(x))\ > - ^)\Y\ > ||y| 

and, by Lemma A. 1, 

Vol { x e R I \x\ < D, \e(X + x,Y + u(x))\ < e } < 4^ < 16 f- 

Case 2: \Y\ < k\X\. In this case we bound the uf l re-derivative away from zero. 
We claim that the dominant term comes from one derivative acting on F and is 2 — 1 
derivatives acting on G. Observe that for \X\, \Y\, \x\ < D with D sufficiently 
small 

\^u(x)\ < 0(D) for0<m<z/ 2 

since u(x) = x U2 c(x, u(x)) and consequently 

\F(x,X,Y)\ < 0(D) 
\^F(x,X,Y)\ > 1- 0(D) 
\^F(x,X,Y)\ < 0(D) forl<m<z/ 2 

Furthermore, since c(X + x, Y + u(x)) — c(x,u(x)) vanishes to first order in 
(X,F),forall x, 

\<§^[x» 2 {c(X + x,Y + u(x)) - c(x,u(x))}]\ < 0(D)[\X\ + \Y\] 

for < m < v 2 and 

\<§^[x» 2 {c(X + x, Y + u(x)) - c(x, u(x))}]\ < 0(1)[\X\ + \Y\] 

for m = v 2 . Since (X + x) 1 " 2 — v 2 Xx V2 ~ l — x 1 " 2 vanishes to second order in X, 
for all x, 

\^[{(X + xy 2 -v 2 Xx» 2 - 1 -x» 2 }c(X + x,Y + u(x))]\ < 0(\X\ 2 ) 

< 0(D)\X\ 
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for all < m < v 2 . Finally 

{^[^Xx^ciX + x,Y + u(x))]\ < 0{D)\X\ for0<m<z/ 2 -l 

\^[v 2 Xx^- l c{X + x,Y + u{x))}\ > k'\X\-0{D)\X\ fox m = v 2 - I 

^[^Xx^ciX + x^Y + u{x))}\ < 0(1)\X\ foxm = v 2 

with k' = v 2 \ inf \c(x, y)\ > 0. Consequently, 

|jg_ e (X + x,Y + u(x))\ 

I dF d v i- x G , (v 2 \ d m F d^- m G | 

\ 2 dx 2.^1 \m) dx m dx v 2~™ I 

m=0 

> (1 - - 0{D))\X\ - 0{D)[\X\ + \Y\] 

> (1 - 0(D))(k - 0(D))\X\ - 0(D)(1 + k)\X\ 



if D is small enough. Hence, by Lemma A.l, 

Vo1{xgM| |x| <D, \e(X + x,Y + u(x))\ <e} 

< o^2+l / e \ 1/ v 2 <- 0^2+1 C 2x/l+K 2 e \ l/"2 



3.2.2 Length of overlap - general case 

Proposition 3.3 Assume H1-H6 with r > 2n + 1. There is a constant D > 
such that for all < 5 < 1 and eac/z sign ± the measure of the set o/p G I 2 such 
that 

!^{ke JP| |e(p ± k)| < ikF }) > (^i) 1 /^' /or some j < 

w most D5 2 . Here £ is the Euclidean measure (length) on T. Recall that no is 
the largest nonflatness or nonnesting order plus one. 

Lemma 3.4 Let r > 2n + 1. For each p6i 2 and k G JF, there are constants 
d,D'>0 (possibly depending on p and k) such that for each sign ±, all j < 0, 
and all pGl 2 obeying |p — p| < d 

*({ke.F||k-k|<d, |e(p±k)|<MO) <D'(^) 1/no 
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Proof of Proposition 3.3, assuming Lemma 3.4. For each pel 2 and k G T, let 
df> t, D'„ r be the constants of the Lemma and set 

Pi K? p,k 

<\k = { (P, k) G M 2 x T | |p - p| < d^, |k - k| < d~ t } 

Since T = { k | e(k) = } and { k | |e(k) | < 1 } are compact, there is an R > 
such that if |p| > R, then {kef |e(p ± k) | < M j } is empty for all j < 0. 
Since { p G IR 2 | |p| < R } x jFis compact, there are (pi, ki), . . ., (p N , k N ) such 
that 

N 

{ P GM 2 ||p|< J R}x^c|J^p l ,k l 

i=i 

Fix any < 5 < 1 and set, for each 1 < i < N, 5 t = (ND'_ k ) n °5. If |p - pi\ > 
Si for all 1 < % < N, then for all j < 

£({kG^| |e(p±k)| <AP'}) 

U {kG^||k- k~^| < d~^, |e(p ± k)| < M J }) 

l<i<n 
|p-Pil<d- . j- 

— ^Pi,ki (|p-pi|) < ^Pi,ki ( <5, ) 



i=l i=l 

AT 

< 



Ei(¥)' / '™ = (¥) 1/ "° 

i=i 

Consequently the measure of the set of p G I 2 for which 

t({ k G T | |e(p ± k)| < M J }) > (If) 1 ^ for some j < 
is at most 

N N 



< where D = ^( ND 'f, i ,^ 



\2n 

n \ i y is - r 

v p l5 ki- 

i=l i=l 



Proof of Lemma 3.4. We give the proof for p + k. The other case is similar. 
In the event that p + k ^ T, there is a d > and an integer j < such that 
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{ k G T | |k - k| < d, |e(p + k)| < M j } is empty for all |p - p| < d and 
j < jo- So we may assume that p + k G JF. 

Case 1: p + k is not a singular point. By a rotation and translation of the k 
plane, we may assume that p + k = and that the tangent line to T at p + k is 
k 2 = 0. Then, as in Section 3.1, there are 

o v G N with 2 < v < n and 

o C r ~ v functions a(q) and 6(q) that are bounded and bounded away from 
zero 

such that 

e(q)=a(q)(fc-?T&(q)) 

in a neighbourhood of 0. (Choose q 2 a(q) = e(q) — e(q)| 92=0 and qia(q)b(q} = 

e(q)| fe =0 .) If the tangent line to T at k (when k is a singular point the tangent 
line of one branch) is not parallel to k\ — and if d is small enough, we can write 
the equation of T for k within a distance d of k (when k is a singular point, the 
equation of the branch under consideration) as k 2 — k 2 = (ki — ki) v ' c(k\ — hi) 
for some 1 < v' < n and some (7 r - n o-i f unc ti n c that is bounded and bounded 
away from zero (when k is a singular point, by Lemma 3.1). Then, for k G T, we 
have, writing p — p = (X, Y) and fci — k\ — x, 

e(p + k) = e(p-p + k-k) = e((X,Y) + (x,x u 'c(x))) 

= A(x, X, Y){Y + x v 'c{x) - (X + x) u B(x, X, Y)) 

where 

A{x,X,Y) = a(X + x,Y + x u 'c(x)) 
B(x,X,Y) = b(X + x,Y + x u 'c(x)) 

and c(x) are bounded and bounded away from zero. 

Observe that y = x v ' c(x) is the equation of a fragment of T translated so as 
to move k to and y = x"b(x, y) is the equation of a fragment of T translated 
so as to move p + k to 0. If p = these two fragments may be identical. That 
is x v c(x) = x u b(x,x u c(x)). (Of course, in this case v = v' .) If p ^ 0, the 
nonnesting condition says that there is an n G N such that if y = x u b(x, y) is 
rewritten in the form y = x v C(x), then the n th derivative of x v C(x) — x v 'c(x) 
must not vanish at x = 0. Let n < n be the smallest such natural number. Since 
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derivatives of x"C(x) at x = of order strictly lower than n agree with the cor- 
responding derivatives of x u 'c(x), the n th derivatives at x = of x u b(x, x v C(x)) 
and x u b(x, x v c(x)) coincide. Since x v C(x) = x u b(x, x u C(x)), the n th derivative 
of x u ' c(x) — x u b(x, x u 'c(x)) = x u ' c(x) — x v B(x, 0, 0) must not vanish at x — 0. 

o If z/ < z/, 

^r(Y +x u ' 'c(x)-(X+xY 'B(x,X, Y))=i/\c{x)+O(d)=i/\c(0)+O(d) 
is uniformly bounded away from zero, if d is small enough. 

o If v 1 > V, 

^(Y + x v 'c(x) - (X + x) v B(x, X, Y)) = -v\ B(0, 0, 0) + O(d) 
is uniformly bounded away from zero, if d is small enough. 

o If v' = v and x v ' c(x) ^ x u B(x, 0, 0), then, as the function x u B(x, 0, 0) - 
(X + x) u B(x, X, Y) vanishes for all x if X = Y = 0, 



d n 
dx 71 



d n 
dx n 



Y + x v 'c{x) - {X + ^^(x, x, y) 
y + x u 'c(x) - x u B(x, 0, 0) 

+[x"B(x, 0, 0) - (X + rr) 1 ^, X, y)] 
= ^(/c(i) - x"B(x, 0, 0)) + 0(|X| + \Y |) 

= £,K'c(x) - o, o))|, =0 + 0(d) + o(\x\ + |y I) 

is uniformly bounded away from zero, if d is small enough. 

o If v' = v and x u ' c(x) = x u B(x, 0, 0) and \Y\ < \X\ 

Y + x u 'c(x) - (X + x) u B(x,X,Y) 
= Y -vXx v - x B(x,XX) 

-{(X + x) v - vXx v ~ x - x u }B(x, X, Y) 
-x u {B(x,X,Y) - B(x,0,0)} 



so that 



dx 



^t(Y + x v 'c[x) - (X + xYB(x, X, Y)) 

-v\ X[B(0, 0, 0) + O(d)} + 0(d)0(\X\ + |y I) 



is bounded away from zero by \v\ \XB(0, 0, 0) |, if d is small enough. 
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In all of the above cases, by Lemma A.l, 



£^{ke^\ |k-k| <d, |e(p + k)| <M j }^j < \]l + c\ 2-0+1 (soM£)V-o 

where u is one of v' , v, n or v — 1, the constant c is the inverse of the infimum 
of a(k), the constant c\ is the maximum slope of T within a distance d of k and p 
is either a constant or a constant times X with X at least a constant times |p — p|. 
There are two remaining possibilities. One is that the tangent line to T at k (when 
k is a singular point the tangent line of one branch) is parallel to k\ = 0. This 
case is easy to handle because the two fragments of T are almost perpendicular, 
so that 

t{{ k G T | |k - k| < d, |e(p + k)| < M j < const M j 

The final possibility is 

o If v' = v and x u ' c(x) = x u B(x, 0, 0) and \Y\ > \X\ 

Y + x u 'c(x)-(X + x) u B(x,X,Y) 

= Y - {(X + xf - x u }B(x, X, Y) - x u {B(x, X, Y) - B(x, 0, 0)} 

so that 

\Y + x u 'c(x) - (X + x) u B(x,X,Y)\ > \Y\ - 0(d)0(\X\ + \Y\) > \\Y\ 

if d is small enough. As a result {kGJF| |k — k| < d, |e(p + k) | < M j } 
is empty if \Y\ is larger than some constant times M J . On the other hand, 
if \Y\ is smaller than a constant times M j , then is larger than some 

I I ' |p-p| o 

constant. 

Case 2: p + k is a singular point. By Lemma 3.1, 

e(p + k + Mq) = a(q)( gi - g^6(q))(g 2 - q^c^)) 

where 2 < i/i, z/ 2 < n are integers, M is a constant, nonsingular matrix and a(k), 
6(k) and c(k) are 

C rr-2-max{i^,^} f unct i ons t hat are bounded and bounded away 

from zero. 

Suppose that the tangent line to M _1 JF at k (when k is a singular point, the 
tangent line of one branch) makes an angle of at most 45° with the rr-axis. Other- 
wise exchange the roles of the q 1 and q 2 coordinates. If d is small enough, we can 



20 



write the equation of T for k within a distance d of k (when k is a singular point, 
the equation of the branch under consideration) as 

(M- X (k - k)) 2 = (M^k - k))fv((M- 1 (k - k))J 

for some 1 < v' < no and some C r ~ n °~ 1 function v that is bounded and bounded 
away from zero. Then, writing M _1 (p — p) = (X, Y) and (M _1 (k — k)) 1 = x 
and assuming that k e T, 

e(p + k) = e(p + k + MM -1 (p-p + k-k)) 

= e(p + k + M(X + x,Y + x u 'v(x))) 

= A(x,X,Y)F(x,X,Y)G(x,X,Y) 

where 

A(x,X,Y) = a(X + x,Y + x u 'v(x)) 

F(x,X,Y) = X + x- (Y + x u 'v(x)) Ul B(x,X,Y) 

G(x,X,Y) = Y + x u 'v(x)-(X + xY 2 C(x,X,Y) 

B(x,X,Y) = b(X + x,Y + x u 'v(x)) 

C(x,X,Y) = c(X + x,Y + x u 'v(x)) 

The functions A(x, X, Y), B(x, X, Y), C{x, X, Y) and v{x) are all C^" - 1 and 
bounded and bounded away from zero. 

As in case 1, y — x u 'v(x) is the equation of a fragment of M _1 JF translated 
so as to move k to and y — x V2 c(x, y) is the equation of a fragment of M _1 JF 
translated so as to move p+k to 0. If p = 0, these two fragments may be identical 
in which case x u 'v(x) = x U2 c(x, x v 'v(x)) and u 2 = v' . This case has already been 
dealt with in Lemma 3.2. Otherwise, the nonnesting condition says that there is 
an n € N such that if y = x u ' 2 c(x, y) is rewritten in the form y = x V2 V{x), then 
the n th derivative of x U2 V(x) — x v v(x) must not vanish at x — 0. Let n < no be 
the smallest such natural number. Since x V2 V(x) = x" 2 c(x, x V2 V(x)) and since 
derivatives at of x V2 V(x) of order lower than n agree with the corresponding 
derivatives of x u v(x), the n th derivative of x u v(x)—x u ' 2 c(x, x v v(x)) = x v v(x) — 
x V2 C(x ) 0, 0) must not vanish at x — 0. So the remaining cases are: 

o \fv'<v 2 , then 

^G(x,X,Y) = u'\v(0) + O(d) 

Since 

F{x, X, Y) = O(d) !L. F (x, X, Y) = 1 + O(d) 
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and applying zero to v' — 1 ^-derivatives to G(x, X, Y) gives 0(d), we have 

£^_F(x, X, X, Y) = (v> + 1)! v(0) + O(d) 

uniformly bounded away from zero, if d is small enough. 

O If V 1 > V 2 , 

^G(x, X, Y) = -v 2 \ C(0, 0, 0) + 0(d) 

Again 

F(x,X,y) = 0(d) ^F(x, X, Y) = 1 + 0(d) 
and applying zero to v 2 — 1 ^-derivatives to X, y) gives 0(d), so that 

£g^F(x, X, Y)G(x, X, y) = -(v 2 + 1)! C(0, 0, 0) + 0(d) 

is uniformly bounded away from zero, if d is small enough. 

o If v' — v 2 and x v 'v(x) ^ x U2 C(x, 0, 0) 

^G(x, X, Y) = % [Y + x»'v(x) - x^C(x, 0, 0) 

+ {x U2 C(x, 0, 0) - (X + x) U2 C(x, x, y)] 
= £t(^t;(x) - tf*C(x, 0, 0))L =0 + O(d) + 0(|X| + |y|) 

and applying strictly fewer than n derivatives gives O(d). As in the last two 
cases 

^tF(x, X, Y)G(x, X, Y) = n^(x»'v(x) - x» 2 C(x, 0, 0)) \ X=Q + O(d) 
is uniformly bounded away from zero, if d is small enough. 
The lemma now follows by Lemma A. 1 , as in Case 1 . ■ 
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4 Regularity 

4.1 The gradient of the self-energy 
4.1.1 The second order contribution 

Let 

where the ultraviolet cutoff U (k) is a smooth compactly supported function that 
is identically one for all k with \ik — e(k)| sufficiently small. We consider the 
value 

F(q) = j d 3 k^d 3 k^d 3 k^ 5(k^ + k& -k® -q)C(kW)C(kM)C(k®) 
of the diagram 7(1) 




k& 

The function V is a second order polynomial in the interaction function v. For de- 
tails, as well as the generalization to frequency-dependent interactions, see [11]. 
For the purposes of the present discussion, all we need is a simple regularity as- 
sumption on V. 

Lemma 4.1 Assume H1-H6. IfV(k^,k^, fc (3) , q) is C\ then F(q) is C 1 in the 
spatial coordinates q. 

Proof: Introduce our standard partition of unity of a neighbourhood of the Fermi 
surface [10, §2.1] 

U(k)=J2f(M- 2j ^k -e(k)\ 2 ) 

j<0 

where f(M- 2j \ik - e(k)| 2 ) vanishes unless M j ~ 2 < \ik - e(k)| < M j . We 
have 

C(k) = Y, C^) where C s {k) = tiM^tegMB 

j<0 
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and 



F(q) = J d 3 k^d 3 k^d 3 k^ 5(k^ + fc(2) _ fc(3) _ q ) 

C il (fc( 1 ))C7 ia (A ; ( 2 ))C7,3(fc( 3 ))V 



JlJ2,i3<0 ' 



Route the external momentum g through the line with smallest (i.e. use the 
delta function to evaluate the integral over the fcW corresponding to the smallest 
and apply V q . Rename the remaining integration variables fcW to and p. 
Permute the indices so that ji < j 2 < J3. If the V q acts on V, the estimate is easy. 
For each fixed j 1 , j 2 and j 3 , 

o the volume of the domain of integration is bounded by a constant times 
\ji\M 2jl \j 2 \M 2j2 , by Lemma 2.3 of [1]. (The k and p components con- 
tribute M jl M j2 to this bound.) 

o and the integrand is bounded by const M~^M~^M~^ Z . 
so that 

J2 I d 3 kd 3 p\C n (k)C n (p)C 3 ,(±k±p±q)V (l V\ 

jl<h<33<0 

< const Y \ji\M jl \j 2 \M h M- j3 < const ^ \ji\ Mh \h\ 

jl<?2<J3<0 jl<?2<0 

< const ^ \ji\ 3 M jl 

ji<0 

is uniformly bounded. So assume that the V q acts on Cj 3 (±k±p±q). The terms 
of interest are now of the form 



d 3 kd 3 p V C il (A;)C j2 (p)V q C j 3(±A; ±p±q) 

with ji < j 2 < j 3 < and 

V q C j3 (±A;±p± g ) = ±[ p ^|_ / ( M -2,3|^ _ e( k ) |2) 

. 2M- 2 « e (k)Ve(k) f u , ■ , - 
+ ifco-U) 7 (M-^3|, fc0 -e(k)|- 



(4) 

fc=±fc±p±<? 



Observe that I V q C i3 (±A;±p±g) I < const M~ 2 ^ since |ifc -e(k)| > const AP 3 
and |e(k)| < const' M js on the support of f(M~ 2j:i \ik — e(k)| 2 ). Choose three 
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small constants rj : f},e > such that < e < < r] < f^§ £ and rj < 
2rao+2^+e • ^ ere n ° * s me mte g er m Proposition 3.3. For example, if n = 3, we 
can choose e = \,rj = ^<^ and fj = ^. 

Reduction 1: For any 77 > 0, it suffices to consider ji < j 2 < J3 < (1 — fj)ji. For 
the remaining terms, we simply bound 

d 3 kd 3 p V C n {k)C n ip) V q C j3 (±Jfe ± p ± g) 

const j d 3 kd 3 pf( llk0 M t )l2 )M-^ /( Mz^ )m-^ M~ 2 ^ 



and 



< const \ix\M jl \j 2 \M h M~ 

\ji\M jl \j 2 \M h M~ 2js 



J3>( 1 -'J)J1 

< const \h\ Mh \h\M h M- j3 - {1 -^ 

jl<j2<ji<0 

< const \M \h\ Mh M - * 1- ^' 1 

jl<32<0 

= const Y liil 3 ^ 1 < 00 
h<o 

Reduction 2: For any 77 > 0, it suffices to consider (k, p) with | ±k±p±q— q| > 
M™ 3 for all singular points q. Let H J3 (k) be the characteristic function of the set 

{ k G R 2 | |k - q| > M vj3 for all singular points q } 

If|±k±p±q — q| < M™ 3 for some singular point q, then | Ve(±k ± p ± q) | < 
const M vj3 and we may bound 



d 3 kd 3 p V C jl (k)C j2 (p)V< l C j3 (±k±p± q) (1 - S js (±k ± p ± q)) 

< const J d 3 kd 3 p /( lik °-^ r )M~ jl /( |ipo ~4 p)|2 )M' j2 M^ 2 '^ 3 

< const|j!|M jl \j 2 \M h M-V-^ 3 

and 
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\ji\ Mh \h\M h M- (2 -^ 3 

ii<i2<i3<o 

< const \ji\M jl \j 2 \M h M-^-^i 

h<j2<o 

< const b'i| 2Mil M- j2+vj2 

ii<i 2 <o 

< const \ji\ 2 M jl M~ jl+Vjl 

h<o 

= const ^ \ji\ 2MVjl < oo 

ji<0 



Current status: It remains to bound 



E | / d 3 kd 3 pVC jl (k)C j MV< l C j3 (±k±p±q)Z j3 (±k±p±q) 

31<32<h<0 
33<(.1-V)il 

< const £ /" rf 3 ^ 3 P /(fcg^)M^ /(J2^1!)M-* 



il<J2<i3 <0 
33<( 1 -'7)J1 



M-^ S j3 (±k ± p ± q)x, 3 (±k ± P ± q) 
< const £ M" 2 ^ I d 2 kd 2 p Xjl (k) Xj2 (p) ( Xj3 S j3 ) (±k ± p ± q) 



Jl<J2<i3<° 



< const Y, M ~ 2n i d2kd2p ^ X » (^Sm) (±k ± p ± q) 

Jl.J2.i'3<° 
<J1 .32 <33 

where Xj (k) is the characteristic function of the set of k's with |e(k) | < M° . 

Make a change of variables with ±k ± p ± q becoming the new k integration 
variable. This gives 

const £ M~ 2js I d 2 M 2 V xn (±k ± p ± q) Xn (p) X* (k) S ja (k) 

ov. ^ 



Jl>J2iJ3<° 
3_ 



A<il,32<33 
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Reduction 3: It suffices to show that, for each fixed k, p and q in IR 2 , there are 
(possibly k, p, q dependent, but ji independent) constants c and C such that 

/ d'k f rf 2 p X , 1 (±k±p± q ) Xj2 (p) X , 3 (k)S J 3(k) 

J\k-k\<c J\p-p\<c 

< C^W^M^-^M^- 11 ''^ (5) 

for all q obeying |q — q| < c. The constant v( will be chosen later and will obey 
e — rj — rj' > 0. Since 

{ (k, p, q) G I 6 |e(k)| < 1, |e(p)| < 1, |e(±k ± p ± q)| < 1 } 

is compact, once this is proven we will have the bound 

const M ~ 2n I d2kd2 P Xii (±k ± P ± q) Xn (p) Xj 3 (k) S is (k) 

< const M- 2h \j 2 \M j2 M^"^ 3 M £h - r '' h 

< const {^-fM^-^ 

33<0 

< const 

since e > rj + rj'. Furthermore, if k or p or ±k ± p ± q does not lie on T, we can 
choose c sufficiently small that the integral of (5) vanishes whenever |q — q| < c 
and \ji | , | j 2 1 , | jz | are large enough. So it suffices to require that k, p and ±k±p±q 
all lie on T. 

Reduction 4: If k is not a singular point, make a change of variables to p = e(k) 
and an "angular" variable 9. So the condition Xj 3 (k) ^ forces \p\ < const M- 73 . 
If k is a singular point, the condition S J3 (k) ^ forces |k — k| > const M™ 3 and 
this, in conjunction with the condition that Xj 3 (k) ^ 0, forces k to lie fairly near 
one of the two branches of T at k at least a distance const M™ 3 from k. Using 
Lemma 3.1, we can make a change of variables such that e(k(p, 9)) = p9 and 
either \9\ > const M™ 3 or \p\ > const M^ 3 . Possibly exchanging the roles of p 
and 9, we may, without loss of generality assume the former. Then the condition 
Xj 3 (k) 7^ forces \p\ < const M^' n ^ 3 . Thus, regardless of whether k is singular 
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or not, 



/ d 2 k / d 2 p Xh (±k ± p ± q) Xh (p) Xh ( k ) E J3 ( k ) 

-/|k-k|<c J\p-p\<c 

- const / | fl |<i dpde I d2 PXjA±Hp,0) ±p±q) xj 2 (p) 



< const / dpd0 / rf 2 p Xj'(± k (0, 0) ± p ± q) x J2 (p) 

^iPi^constMd-^s ■ / Ip-pI<c 

< const M* 1 -")'' 3 / d0 / d 2 p X /(±k(O,0)±p±q)x J2 (p) 
V|6>I<1 ^Ip-pKc 



'I"I<1 ^|P-P|<C 

where AP" = AP' 1 + const M^-^' 3 < const M^ 1-7 ^ 3 . Thus it suffices to prove 
that 

/ d 2 p [ d9 X A±H0,0)±p±q)XjM<C\j2\M^M^'^ 

^|p-p|<c ^|6>|<1 

for all q obeying |q — q| < c. 

The Final Step: We apply Proposition 3.3 with j = f and 5 = M^+^Z 2 . 
Denote by x(p) the characteristic function of the set of p's with 



//({ -l<0<l||e(±p±q±k(O,0))|<M''}) >*(mV 



Mi \ 1/np 
M( 1+e: W3/ 2 / 



where c\ is the supremum of j- (s is arc length). If ±p ± q is not in the set of 
measure D5 2 specified in Proposition 3.3, then 



l/no 



*({ k G T | |e( ± p ± q ± k)) | < M j ' }) < (*f ) 

and hence 



/i({ - 1 < < 1 | |e(±p±q±k(O,0))| < M j ' } 

= J d6 X (\e{ ± p ± q ± k(0, 9))\ < ilP') 



= J dsgx(|e(±P±q±k)|<M'"') 
<Ci^({ke T\ |e(±p±q±k))| < M j ' }) 



< ri f ^/"o 

Cl V M (l+ £ )j 3 /2 J 
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so that x(p) = 0. Thus x(p) vanishes except on a set of measure DM (l+£ ^ 3 and 
/ d 2 p [ dQ x,'(±k(0, 9) ± p ± q) X J2 (P) 

J ./|0|<1 

< / d 2 p x(p)XjM [ d9 X/ (±k(0, 0) ± p ± q) 
J J\e\<i 

+ / d 2 p (1 - X(P))X J2 (P) / d0x,v(±k(O,0)±p±q) 

J -/|6>|<1 

<2 | d 2 p x(p) + const | ^PX^pX^^,^) 17 " 

< const M( 1+£ ^+ const | J ' 2 |M^( ]gil ^) 1/no 

< const M j3 M £j:i + const |j 2 |M J2 (M (1 -^ 3 - i ^ 3 ) 1/no 

< const M h M {£ ~^ )j3 + const |j 2 1 M J2 (M ( 5"'?- f)«) 1/no 

since j 3 < (1 - fj)j 2 < j 2 - j^js 

< const \j 2 \M j2 M £j3 - v ' j3 

provided e < > and r/ > ^ + 

We choose ?/ = — + Then the remaining conditions 

^ + V s>*7 + £ + ^ 2 S^> a fl ??<g^ ( 6 ) 

and 

„' > J!_ f, < jL_ = 2r i+ £ 

'I — 1-77 ' — !+»/ 2n +27j+£ 

are satisfied because of the conditions we imposed when we chose r/, fj and e just 
before Reduction 1 . 

We have now verified that V q F(g) is a uniformly convergent sum (over j u j 2 , 
and j 3 ) of the continuous functions V q / d 3 kd 3 p V C jl (k)C j2 {p)Cj. A (±k ±p± q). 
Hence F(q) is C 1 in q. ■ 



4.1.2 The general diagram 

The argument of the last section applies equally well to general diagrams. 

Lemma 4.2 Let G(q) be the value of any two-legged 1PI graph with external 
momentum q. Then G(q) is C 1 with respect to the spatial components q. 
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Proof: We shall simply merge the argument of the last section with the general 
bounding argument of [1, Appendix A]. This is a good time to read that Appendix, 
since we shall just explain the modifications to be made to it. In addition to the 
small constants 77, 7/, 77, e > of Lemma 4.1, we choose a small constant e > 
and require 1 that 

°< £ Z£;,0<V<%&e, V ' = % + ^and^<mm{ V ',e- V - V '} 
and 

e < min{r7 , 77 , (1 + e — rj — rj')(l — fj) — 1} 

with n being the integer in Proposition 3.3. All of these conditions may be satis- 
fied by 

o choosing < e < and then 

o choosing < 77 < ^°~j, £ (by (6), this ensures that e — 77 — rj' > 0) and then 

o choosing < fj < 1 so that < min {77', e — 77 — 77'} (this ensures that 
the expression (1 + e — 77 — 77') (1 — fj) — 1 > 0) and then 

o choosing e > so that e < min{77 , 77 , (1 + e — 77 — 77') (1 — fj) — 1} 

As in [1, Appendix A], use [1, (22)] to introduce a scale expansion for each 
propagator and express G(q) in terms of a renormalized tree expansion [1, (24)]. 
We shall prove, by induction on the depth, D, of G J , the bound 

sup|^«^ 1 G' J (g)| < const n |jf n - 2 M J M~ SoJ M" Sl(1 - £ " )j (7) 
JeJ(j,t,R,G) q 

for s , si G {0, 1}. The notation is as in [1, Appendix A]: n is the number of 
vertices in G and J(j, t, R, G) is the set of all assignments J of scales to the 
lines of G that have root scale j, that give forest t and that are compatible with 
the assignment R of renormalization labels to the two-legged forks of t. (This is 
explained in more detail just before [1, (24)].) If s — and si = 1, the right hand 
side becomes constJj| 3n-2 M £:? , which is summable over j < 0, implying that 
G(q) is C 1 with respect to the spatial components q. If s 1 = 0, (7) is contained in 
[1, Proposition A. 1], so it suffices to consider si — 1. 

'The first three conditions as well as the condition that < rj were already present in 
Lemma 4. 1 . The other conditions are new. 
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As in [1, Appendix A], if D > 0, decompose the tree t into a pruned tree t and 
insertion subtrees r 1 , • • • , r m by cutting the branches beneath all minimal Ef = 2 
forks /i, • • • , f m . In other words each of the forks /i, • • • , f m is an Ef = 2 fork 
having no Ef = 2 forks, except 0, below it in t. Each consists of the fork fa and 
all of t that is above /j. It has depth at most — 1 so the corresponding subgraph 
Gf i obeys (7). Think of each subgraph Gf i as a generalized vertex in the graph 
G = G/jG/j, • • • , Gf m }. Thus G now has two as well as four-legged vertices. 
These two-legged vertices have kernels of the form Tj(/c) = V <■ IGf^k) 
when fi is a c-fork and of the form Ti(k) = ^2- > (11 — tjGf^k) when /j is 
an r-fork. At least one of the external lines of Gf i must be of scale precisely j 7r (/ l ) 
so the momentum k passing through Gf z lies in the support of Cj. f , y In the case 
of a c-fork / = f { we have, as in [1, (27)] and using the same notation, by the 
inductive hypothesis, 



E E sup d^£G J f f (k) < J2const nf \jf\ 3n f- 2 M^M~ S1 ^- 



£ )if 



jf<Mf) J f eJ(j f ,t f ,R f ,G f ) k ' ' jf<Mf) 



< const^l^^l^-^^^M- 8 ^ 1 -^'^/) (8) 



for s 1 — 0, 1. As lG J f (k) is independent of k derivatives with respect to k Q may 
not act on it. In the case of an r-fork / = f h we have, as in [1, (29)], 



E E su p + °) - tyG J /(k) 



- E E M^ 1 -' ^/) sup d ko d^G/(k) 

jf>Mf) J f eJUf,tf,Rf,G f ) 

< const n/ M (1 - so) ^(/) E Ij/I 3 "'"^- 81 ^-^'' 

< a)nst n J^(/)|^ / " 1 M^(«M- ao ^(«M -ai(1 -^(« (9) 

Denote by J the restriction to G of the scale assignment J. We bound G J , 
which again is of the form [1, (31)], by a variant of the six step procedure followed 
in [1, Appendix A]. In fact the first four steps are identical. 

1 . Choose a spanning tree T for G with the property that Tfl Gj is a connected 
tree for every / e t(G J ). 
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2. Apply any g-derivatives. By the product rule each derivative may act on 
any line or vertex on the "external momentum path". It suffices to consider 
any one such action. 

3. Bound each two-legged renormalized subgraph (i.e. r-fork) by (9) and 
each two-legged counterterm (i.e. c-fork) by (8). Observe that when s' Q 
/co-derivatives and s[ k-derivatives act on the vertex, the bound is no worse 
than M~ s 'o j M~ s i ( 1_£ W times the bound with no derivatives, because we nec- 
essarily have j < jn(f) < 0. 

4. Bound all remaining vertex functions, u Y , (suitably differentiated) by their 
suprema in momentum space. 

We have already observed that if si = 0, the bound (7) is contained in [1, Propo- 
sition A.l], with s = 0. In the event that si = 1, but the spatial gradient acts on a 
vertex, [1, Proposition A.l], again with s = but with either one v replaced by its 
gradient or with an extra factor of M~^~ 6 ^ coming from Step 3, again gives (7). 

So it suffices to consider the case that si = 1 and the spatial gradient acts on 
a propagator of the "external momentum path". It is in this case that we apply the 
arguments of Lemma 4.2. The heart of those arguments was the observation that, 
when the gradient acted on a line £ 3 of scale j 3 , the line £ 3 also lay on distinct 
momentum loops, and A^ 2 , generated by lines, £ 1 and £ 2 of scales j £l and jg 2 
with jg 1 , j( 2 < je 3 . This is still the case and is proven in Lemmas 4.3 and 4.4 
below. So we may now apply the procedure of Lemma 4.1. 

Reduction 1: It suffices to consider j < je s < (1 — fj)j. For the remaining terms, 
we simply bound the differentiated propagator, as in the argument following (4), 
by 

\d%dqC jta (k t3 (q))\ < const M- 2 ^M- s o% 

< const M- S '" je s M- (1 -^' 

< const M~ s 'o je 3 M~ {1 - £ ~ )j 

So for the terms with > (1 — f})j, the effect of the spatial gradient is to degrade 
the si — bound by at most a factor of const M~^ l ~ 6 ^ and we may apply the 
rest of [1, Proposition A.l], starting with step 5, without further modification. 
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Reduction 2: It suffices to consider loop momenta (ki) ee Q,f, in the domain of 
integration of [1, (31)], for which the momentum k^ 3 flowing through £ 3 (which 
is a linear combination of q and various loop momenta) remains a distance at 
least M^3 from all singular points. If k^ 3 is at most a distance M v ^3 from some 
singular point then |<9ke(k^ 3 )| < const M vji 3 and we may use the Ve in the 
numerator of (4) to improve the bound on the £ 3 propagator to 



Once again, in this case, the effect of the spatial gradient is to degrade the si = 
bound by at most a factor of const M~^~ £ ^ j and we may apply the rest of the 
proof [1, Proposition A.l], starting with step 5, without further modification. 

Reduction 3: Now apply step 5, that is, bound every propagator. The extra <9 q 
acting on Cj e gives a factor of M~ sije 3 < M~ Slj worse than the bound achieved 
in step 5 of [1, Proposition A.l]. Prepare for the application of step 6, the inte- 
gration over loop momenta, by ordering the integrals in such a way that the two 
integrals executed first (that is the two innermost integrals) are those over ke 1 and 
ki 2 . The momentum flowing through £ 3 is of the form ke a = ±k £l ± k i2 + q', 
where q' is some linear combination of the external momentum q and possibly 
various other loop momenta. Make a change of variables from and k^ 2 to 
k = k^ 3 = ±k^ ± k^ 2 + q' and p = k^ 2 . It now suffices to show that, for each 
fixed k, p and q' in IR 2 , there are (possibly k, p, q' dependent, but j e . independent) 
constants c and C such that 



for all q' obeying |q' — q'| < c. Recall that Xj(k) is the characteristic function of 
the set of k's with |e(k) | < M j and Sj(k) is the characteristic function of the set 



Once proven, the bound (10) (together with the usual compactness argument) 
replaces the bound 



d qo d <i C h 3 (**s (?) ) I < const M ~ 2jt3 M~ s >3 M™3 

< const M- je 3 M- S '" je 3 M- {1 - T,)j 

< const M~ je 3 M- s >3 M~ (1 - £ ~ )j 




< C \je 2 1 M je * M (1 "*'3 M £k 3 (10) 



{ k G R 2 | |k - q| > M vj for all singular points q } 




d 2 k £l / d 2 k h Xje AK) X je AK) < const \j tl \M»i \j t2 \M* 
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used in step 6 of [1, Proposition A.l]. Since j < je 1 ,je 2 < 3e 3 < (1 — the 
bound (10) constitutes an improvement by a factor of 

const M je i M (1 ~ v)je 3 M £je 3 M~ jt i M~ k * < const M^^M"^ 

< const mK 1 -"-"'-^)^-^- 1 ^' 

< const M £_i 

So once we proven the bound (10), we may continue with the rest of the proof 
of [1, Proposition A.l], without further modification, and show that the inductive 
hypothesis (7) is indeed preserved. 

In fact the bound (10) has already been proven in Reduction 4 and The Final 
Step of the proof of Lemma 4.1. ■ 

Lemma 4.3 Let G be any two-legged 1 PI graph with each vertex having an even 
number of legs. Let T be a spanning tree for G. Assume that the two external 
legs of G are hooked to two distinct vertices and that £3 is a line of G that is in 
the linear subtree of T joining the external legs. Recall that any line £ not in T 
is associated with a loop Ae that consists of t and the linear subtree ofT joining 
the vertices at the ends of £. There exist two lines £\ and £ 2 , not in T such that 

£ 3 e A h n A h . 

Proof: Since T is a tree, T\ {£3} necessarily contains exactly two connected 
components Ti and T 2 (though one could consist of just a single vertex). On the 
other hand, since G is 1PI, G \ {£3} must be connected. So there must be a path 
in G \ T that connects the two components of T \ {^3}. Since T is a spanning 
tree, every line of G \ T joins two vertices of T, so we may alway choose the 
connecting path to consist of a single line. Let £ x be any such line. Then £3 e . 
If G \ {£1, £3} is still connected, then there must be a second line £ 2 7^ £\ of G \ T 
that connects the two components of T \ {^3}. Again £ 2 £ A ei . If G \ {£1, £3} 
is not connected, it consists of two connected components G\ and G 2 with G\ 
containing 7\ and G 2 containing T 2 . Each of Ti and T 2 must contain exactly one 
external vertex of G. So each of G± and G 2 must have exactly one external leg 
that is also an external leg of G. As £\ and £3 are the remaining external legs of 
both G\ and G 2 , each has three external legs, which is impossible. ■ 

Lemma 4.4 Let G and T be as in Lemma 4.3. Let J be an assignment of scales 
to the lines of G such that T n Gj is a connected tree for every fork f G t{G J ). 
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(See step 1 of [1, Proposition A. 1 ].) Let £ 3 G T. Let £ G G \T connect the two 
components ofT \ {£ 3 }, then j e < jg 3 . 

Proof: Let G' be the connected component containing £ of the subgraph of G 
consisting of lines having scales j > j t . By hypothesis, G' n T is a spanning tree 
for G' . So the linear subtree of T joining the vertices of £ is completely contained 
in G'. But £ 3 is a member of that line's subtree and so has scale je 3 > je. ■ 

4.2 The frequency derivative of the self-energy 

We show to all orders that singularities in this derivative can only occur at the 
closure of the lattice generated by the Van Hove points. That the singularities 
really occur is shown by explicit calculations in model cases in the following 
sections. 

Lemma 4.5 Let G(q) be the value of any two-legged 1PI graph with external 
momentum q. Let B be the closure of the set of momenta of the form (0, q) with 
q = Yli=i(~ where n G N and, for each 1 < i < n, Sj G {0, 1} and q« is a 
singular point. Then G(q) is C 1 with respect to q on M, 3 \ B. 

Proof: The proof is similar to that of Lemma 4.2. Introduce scales in the 
standard way and denote the root scale j. Choose a spanning tree T in the standard 
way. View c- and r-forks as vertices. The external momentum is always routed 
through the spanning tree so the derivative may only act on vertices and on lines of 
the spanning tree. The cases in which the derivative acts on an interaction vertex 
or c-fork are trivial. If the derivative acts on an r-fork, the effect on the bound 
[1, (29)], namely a factor of M~ j <f>, is the same as the effect on the bound [1, 
(32)] when the derivative acts on a propagator attached to the r-fork. So suppose 
that the derivative acts on a line £ 3 of the spanning tree that has scale j 3 . We know 
from Lemma 4.3 in the last section that there exist two different lines £\ and £ 2 , 
not in T such that £ 3 lies on the loops associated to £\ and £ 2 . We also know, from 
Lemma 4.4, that the scales j' of all loop momenta running through £ 3 , including 
the scales ji and j 2 of the two lines chosen, obey j' < j 3 . 

Reduction 1: It suffices to consider j < ji,j 2 < J3 < (l — fj)j. For the remaining 
terms, we simply bound 

I j-C j3 (±q ± internal momenta) | < const M~ 2h < const M~ h M^ 1 ^ 
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After one sums over all scales except the root scale j, one ends up with 

const \j\ n M j M-^-^ j = const \j\ n M^ j 
which is still summable. 

Reduction 2: Denote by k, p and ±k ± p + q' the momenta flowing in the lines 
l\, £ 2 and £3, respectively. Here q' is plus or minus the external momentum q 
possibly plus or minus some some other loop momenta. In this reduction we 
prove that it suffices to consider (k, p) with |k — k| < and |p — p| < 
and I ± k ± p ± q' — q| < for some singular points k, p and q. 

Suppose that at least one of k, p, ±k ± p ± q' is farther than M r ' j from 
all singular points. We can make a change of variables (just for the purposes 
of computing the volume of the domain of integration) such that k is at least a 
distance from all singular points. After the change of variables, the indices 
ji, J2 and J3 are no longer ordered, but all are still between j and (1 — fj)j. Let 
Hj (k) be the characteristic function of the set 

{keK 2 ||k-k|> M vj for all singular points k } 

We claim that 

Vol{(k,p) G R 4 |x jl (k)S J (k) Xj2 (p) Xj3 (±k±p±q')^0} 

< C\j 2 \M h ^ j M j2 M £j -^ 

This would constitute a volume improvement of M^ -2 ' 7 ^ and would provide 
summability if 2rj < e. By the usual compactness arguments, it suffices to show 
that, for each fixed k, p and q in R 2 , there are (possibly k, p, q dependent, but 
j, ji independent) constants c and C such that 

/ d 2 k [ rf 2 p Xjl (k)S,(k) Xj2 (p) X , 3 (±k±p±qO<C|j 2 |M^-^M^M^ 

^|k-k|<c ^|p-p|<c 

for all q' obeying |q' — q| < c. Furthermore, if k or p or ±k ± p ± q' does not 
lie on T, we can choose c sufficiently small that the integral vanishes whenever 
|q' — q| < c and |ji|, |j 2 |, |j 3 | are large enough. So it suffices to require that k, p 
and ±k ± p ± q' all lie on T . 

If k is not a singular point, make a change of variable to p = e(k) and an 
"angular" variable 9. If k is a singular point, the condition 2j(k) 7^ forces 
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|k — k| > M V K We can make a change of variables such that e(k(p, 6)) = p9 and 
either \9\ > const or \p\ > const M™. Possibly exchanging the roles of p 
and 6, we may, without loss of generality assume the former. In all of these cases, 
the condition Xji(k) 7^ forces \p\ < const M jl ~ vj . Thus 

f d 2 k [ rf 2 p X , 1 (k)S J (k) X , 2 (p) X , 3 (±k±p±q') 

J\k-k\<c J\p-p\<c 

- COnSt / ,,,<! dpd0 / rf2p Xj 2 (p)x J3 (±Hp, 0) ± p ± q') 
dpdfl / rf 2 px / (±k(0^)±p±q')x J2 (p) 

|p|<con st M«-™ J \P-P\<C 



< const 



< const M jl ~ vj f d6 f d 2 p Xj '(±k(0, 0) ± p ± q') Xj2 (p) 
J|e|<i ^|p-p|<c 

where M j ' = M h + const M jl ~^ < const M^'^ . Thus it suffices to prove 
that 

/ rf 2 p / ^ X,-'(±k(0, 6)±p± q') X J2 (P) < C\j 2 \M*M*-* 
J\v-v\<c J\e\<i 



for all q' obeying |q' — q| < c. 

We again apply Proposition 3.3 with j = f and 5 = M^ 1+e ^ 2 / 2 . If we denote 
by x(p) the characteristic function of the set of p's with 



A*({ -1<0<1\ |e(±p±q'±k(O,0))| < M j ' }) > ci( 



M'' \l/np 

M (l+ £ )j 2 /2 J 



where c x is the supremum of ^ (s is arc length), then x(p) vanishes except on a 
set of measure DM^ l+£ ^ and 

/ d 2 P [ de xA±p ± q' ± k(o, e)) Xj2 (p) 

J ./|0|<1 

< / ^ 2 p x(p)x J2 (p) / ^ xr(±p ± q ± k(o, ^)) 

J J\e\<i 

+ [ d 2 p (1 - x(p))XniP) [ d9 X f(±P±4±k(0,9)) 
J J\o\<i 

<2j d 2 p x(p) + const | d 2 px,(p)(^) 1/r, ° 
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< const const ^M^i^M^) 1 ^ 

< const M j2 M £{1 - f ' )j + const |j2|M J ' 2 (M( 1 -"-^'- 1 * £j ') 1/no 

< const M j2 M £ ^~^ j + const \j 2 \M h M^~ v)j 

< const |j 2 |M J ' a M ej ' - '» 

provided fj < rj < ^ and £ < min {|, 

End stage of proof The momentum flowing through the differentiated line is of 
the form ±g± k\ ± . . . ± k n where q is the external momentum and each of the fcj's 
is a loop momentum of a line not in the spanning tree whose scale is no closer to 
zero than the scale of the differentiated line. If q £ B, then either q is nonzero or 
the infimum of | ±q±ki ±. . .±k„ — k n+ i], with the kj's running over all singular 
points, is nonzero. So there exists a j such that when the root scale obeys j < j , 
either the zero component of one of k\, . . ., k n , ±q ± k± ± . . . ± k n has magnitude 
larger than const M\ in which case the corresponding covariance vanishes, or 
the distance of one of k l5 . . ., k„, ±q ± ki ± . . . ± k ra to the nearest singular point 
is at least M r " , in which case we can apply the argument of reduction 2. (If it is 
one of ki, . . ., k n whose distance to the nearest singular point is at least M vj , we 
may choose as the l\ of Lemma 4.3 the line initiating that loop momentum.) ■ 



5 Singularities 

In this section, we do a two-loop calculation for a typical case to show that the 
frequency derivative of the self-energy is indeed divergent in typical situations, 
and to calculate the second spatial derivative at the singular points. 

By the Morse lemma, there are coordinates (x, y) such that in a neighbourhood 
of the Van Hove singularity, the dispersion relation becomes 

e(k) = e(x,y) =xy. 

Here we consider the case of a Van Hove singularity at k = with e(k) = 
(In particular the Van Hove singularity is on the Fermi surface for \i = 0.) The 
nonlinearities induced by the changes of variables are absent in this example, and 
moreover, the curvature is zero on the Fermi surface. We rescale to a; = I 1 , 
y = I 2 - and, for definiteness, take the integration region for each variable to be 
[-1, if 
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For this case, we determine the asymptotics of derivatives of the two-loop con- 
tribution to the self-energy as a function of q for small q . We find that again, the 
gradient of the self-energy is bounded (in fact, the correction is zero in that case) 
but that the g -derivative is indeed divergent. Power counting by standard scales 
suggests that this derivative diverges at zero temperature like | log go I • However, 
there is a cancellation of the leading singularity which is not seen when taking ab- 
solute values, so that the true behaviour is only (logg ) 2 . We then also calculate 
the asymptotics of the second spatial derivative and find that it is of the same order 
as the first frequency derivative. Finally, we do the calculation for the one-loop 
contributions to the four-point function, to compare the coefficients of different 
divergences in perturbation theory. 

The physical significance of these results will be discussed in Section 6. 

5.1 Preparations 

We restrict to a local potential. Since we have only considered short-range inter- 
actions, the potential is smooth in momentum space. For differentiability ques- 
tions, a momentum dependence could only make a difference if the potential van- 
ished at the singular points or other special points, so the restriction to a local 
potential, which is constant in momentum space, is not a loss of generality. 

There are two graphs contributing, one of vertex correction type and the other 
of vacuum polarization type (the graphs with insertion of first-order self-energy 
graphs have been eliminated by renormalization through a shift in //). 





The latter gets a (—1) from the fermion loop and a 2 from the spin sum. Thus the 
total contribution is 

S 2(<?o, q) = ~j2 ^2 Wf J d 2 k 1 d 2 k 2 d 2 k 3 S(q - k ± + k 2 - fa) 

C(wi, e(fci)) C(u 2 , e(k 2 )) C(q -uj ± + u 2 , e(fa)) 



with 



C(u,E) 



iu — E 
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Call Ei = e(ki) and 

(F) g = J d Pq (k) F(k) 

where k = (fci, k 2 , k 3 ) and dp q (k) = / d 2 k\d 2 k 2 d 2 k 3 5(q — ki + k 2 — k 3 ). 
The frequency summation gives 

S2(g °' g) -"\ ig + (E 2 -E 3 - ET) ) q (U) 

where fp(E) = (1 + e^)" 1 is the Fermi function and bp(E) = (e pE - l)' 1 is the 
Bose function. Since 

bfi(E 2 - E 3 ) [fp(E 2 ) - f p (E 3 )] = fp{E 2 )[f p {E 3 ) - 1] 

the numerator is bounded in magnitude by 2. The denominator is bounded below 
in magnitude by |g |, so the integrand is C°° in q , for all go 7^ and all (3 > 0. 
Because the integral is over a compact region of momenta, the same holds for 
the integral. The limit (3 — > oo exists and has the same property. The structure 
of the denominator may suggest that for it to almost vanish requires only the 
combination E 2 — E 3 — E x to get small, but a closer look reveals that each E,- L has to 
be small: at T = 0, the Fermi functions become step functions, fp(E) — > Q(—E) 
and bp(E) — > — ©(— E), and then the factors in the numerator imply that all 
summands in E 2 — E 3 — E 1 really have the same sign, i.e. \E 2 — E 3 — Ei\ = 
| E 2 1 + | E 3 1 + | Ei | , so all | Ei \ must be small for the energy difference to be small. 
At finite (3, when the E'j's are "of the wrong sign" the exponential suppression 
provided by the numerator compensates for the |g | > j in the denominator. 
Because 

ig -e(g)-A 2 S 2 (g ,g) = igo(l + iA 2 9 S 2 (0,g))-( e (g)+A 2 S 2 (0,g)) + ... (12) 

and because <9 E 2 (0, q) is purely imaginary, we are interested in 

Z 2 (q) = (l-A 2 Im<9oS 2 (0,g)) _1 . (13) 

The value q = is not an allowed fermionic Matsubara frequency at T > 0. We 
shall keep q ^ in the calculations. In discussions about temperature depen- 
dence, we shall replace |g | by ir/{3. 

The second order contribution to Im <9 S is 



Im<9 £ 2 (g ,g) = {$ qo (E 2 - E 3 - Ej 

IfpiEi) + b p (E 2 - E 3 )] [fp(E 2 ) - fp{E 3 )] ) (14) 
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with 

J ? - ql 

(ig + eY (e* + ql) 



Re 7— 7^ = 7T7^2- d5) 



5.2 go-derivative 

The above expression (14) for Im <9 £ 2 (g , <?) shows that it is an even function of 
q and that q serves as a regulator so that even at zero temperature, a singularity 
can develop only in the limit q — > 0. We therefore calculate it as a function of q 
at zero temperature. In the following, we take q > 0. 

As (3 — > oo, the Fermi function //j(.E') — > Q(—E) and for £ ^ 0, the Bose 
function — > — ©(— -E 1 ). In this limit, the integrand vanishes except when 
E 2 E 3 < and E 1 (E 2 — E 3 ) < 0. This reduces to the two cases 

E 1 > 0, E 2 < 0, E 3 > and E 1 <0,E 2 > 0, £ 3 < 0. (16) 

In both cases, the combination of //3,6/3's in the numerator is —1. Thus 

Im <9 £ 2 (g , q) = -(^90(^2 - E 3 - E ± ) [ 1 > A £ 2 < A £3 > 0) 

+ 1(£ 1 <0AE 2 >0A£ 3 <0) 1 \ 

J / <? 

Recall that we are considering the case of a Van Hove point at k — (x, y) — 
for the dispersion relation e(k) = xy with x = ki/2-n and y = k 2 /2n, so that 
= dx dy. Moreover we set q — 0, and use the delta function to fix E\ in 
terms of -E 2 and E 3 , so that 

£ 2 = xy E 3 = x'y' E 1 = (x - x') (y - y') 

and 

e = E 2 — E 3 — Ei = xy' — Ix'y + x'y. 

Recall also that, for definiteness, we are taking the integration region for each 
variable to be [—1,1]. The sign conditions on the Ei impose conditions on the 
variables x, . . ., which are listed in Appendix B, and which we use to transform 
the integration region to [0, l] 4 . 

At T = 0, only n e M = {1, 2, 3, 4, 9, 10, 11, 12} from the table in Appendix 
B contribute. Thus 

Im<9 £ 2 (g ,0) = - / dxdydx'dy' ^ l{p n ) $ 9o (e„(x,y,x',y')) 



neM 
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with e n and p n given in the table in Appendix B. By (31), (32) and (33), and 
because $ 90 is even in e, all eight terms give the same contribution. Hence 

Im<9 £ 2 (g ,0) = -2I(q ) (17) 

where 

fl pi fl rx' 

/(go) =4/ dy dy' dx' dx $ ?0 ((2x' - x)y' + yx') . (18) 
Jo Jo Jo Jo 

Lemma 5.1 Let I be defined as in (18) and < go < \. Then 

/m9 S 2 (g ,0) = -4 log 2 |logg | 2 - 2Ci |logg | + B(q ) 
where B is a bounded function ofqo, and 

Proof: By (17), it suffices to show that 

J(g ) = 2 log 2 lloggol' + Ci |logg | + B(q ) 

with bounded B . We rewrite the argument e of $ ?0 as e = x'(2y' + y) — xy' . We 
first bound the contribution of y' < q . To do this, bound 



1 

q 2 o~+^ 2 

Use that this bound is decreasing in e and that e > x'y. Therefore 

nl nq pi rx' 

4 dy dy' dx' / dx $ ?0 ((2x - x)y' + yx') 
Jo Jo Jo Jo 

< 4 / dy / dy' dx' / dx 2 

Jo Jo 7o Jo % + \ x y) 

= 4g / dx' / dy 



jo jo ~* ?o + (aty) 2 

„ f 1 1 / 1 

4g dx — arctan — 
Jo Qo Qo 



< 2tt. 
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Thus it suffices to calculate the asymptotic behaviour of 



1 rl 



7(g ) = 4 / dy dy' dx' / dx $ go (x'(2y' + y) - xy') 

JO Jqo JO 



for small g > 0. By (15) 



(e) 



9 £ 



cte g^ + e 2 ' 



(19) 



Because §§ = -y' 



/ dx % {e{x)) = -, 
Jo y 



e 2 + ql 



The integral over x' can now be done, using 



dx' 



ax 



Thus 
with 



(ax 1 ) 2 + ql 2a 

/(go) = h ~ h 



(y'+y)x' 
(2y'+y)x' 



a 



?o 



y' 7o 2{jy' + v) 



log 1 + 



% 2 



2 ^ 
'i »7 



d?7 r^ +jr) d£ 



3 V 



-f log (i + e) 



where we have made the change of variables £ = , d£ = -^dy followed by 
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the change of variables rj — ^, drj — ^dy'. Thus 



Hod) = 



dr) 



i V 



f q o d V / 

Ji 7 ^ 



\VM J [q^ 1 +r 1 ,q- 1 +2r ) \ 



) 



with 



J a =2 



[ ^ log (1 + £ 2 ) > for A C [0, oo) 
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We have 

rb 2 



J M =2 / f log(l + e 2 )= / y log(l + t). 

</a 4 Jo? 1 



(20) 



In our case both integration intervals for J are subsets of [1, oo), so we can expand 
the logarithm, to get 



J[a,b] = f^J [logW+log(l + ^ 

= \ [(logfe 2 ) 2 - (log a 2 ) 2 ] - E ^ ( Q_2n " r ") 

n>l 

= 2 [(log 6) 2 - (log a) 2 ] - £ (a- 2 « - 6" 2 «) 



n>l 



= 2 log(a&) log - - V ( -*» _ 6-2") (21) 

n>l 

The final integral over r\ gives, for the first term, 

f" 1 dn 

/ — J[ V ,2 V ] = 21og2 |loggo| 2 + 2(log2) 2 | log g 1 

-sE^a-oo-,?') 

n>l 

with the last term analytic, and hence bounded, for |g | < I. The second term 
gives two contributions: 

I ^VW+ai] = W-R 

with (here M = q^ 1 ) 
R= *(V), X(rj) = £ t^L [ (M + r?)" 2 " - (M + 217)-* 1 ] 



and 



fM 

W = 2 ' 
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For \q \ < 1 the series for X converges absolutely and gives 

x = E ti n K 1 + ™y 2n - ^ + ^) _an ] 



n>l 

so 1^1 < E„>i go"/^ 2 < 2 9o, hence for |g | < 1, 

R < 2ql\ lng | < go 
In W, scaling back to x = q rj gives 
dx 



W 



2 / — 

x 



a=2 
a=l 



(\ogq 1 + log(l + ax))' 

f 1 dx 2 
= 2/ — [2\ogq 1 \og(l + ax) + (\og(l + ax)) 2 ] a a Z 2 1 . 

Jq X 

Because ax > 0, < log(l + ax) < ax, so W is bounded by a constant times 
l°g 1 - The integral of the same function from to q is of order q \ \ogq \, 
therefore 



[ % log (I±§) + m 



go J 



where 5 is a bounded function. 



5.3 First spatial derivatives 

In our model case, we have E 2 = xy and E 3 = x'y' . Moreover we take ^ = (£, 77) 
so that, fixing k\ by momentum conservation, E 1 = (£ + x — x')(r) + y — y'). It is 
clear from (11) that the spatial derivatives also act on the Fermi function fp(Ei), 
so we get two terms. The derivative of the Fermi function fp(x) is minus the 
approximate delta 6p (x) = (3/ (A cosh 2 ((3x/2)) so that 



-diE 2 (g , g) = ( {die 



r / p J2-/3 . (/l + b 23 )(f 2 — /3J 

-^l-c'ij- : H 



igo + ^ 



(igo + e) 2 



(22) 



where <%e has argument q + k 2 - k 3 , fa = fp(Ei), b 23 = bp(E 2 - E 3 ) and, as 
before, e = E 2 - E 3 - E x . At q = 



9 , 
-^S 2 (g ,0) 



go) +S 2 {I3, go) 
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with 

yl) Mxy) - Mx'y') 
(S,)((x - x')(y - y')) 

y } Cm + e) 2 

[f^ x - x ')(y-y')) + bp(xy-x>y')] 

Here e = xy — x'y' — (x — x')(y — y') = xy' + x'y — 2x'y'. 

Now consider Si and apply the reflection (x, y, x', y') — > (— x, —y, —x', —y') 
to the integration variables. The domain of integration is invariant. The only 
noninvariant factor is y — y' and it changes its sign. Thus Si vanishes. By the 
same argument, S 2 vanishes as well. By symmetry, the same holds for the rj- 
derivative. Thus 

VS 2 (g ,0) = 0. 



Si(P,qo) = I dxdydxdy (y 
-i,i] 4 



S2(P,qo) = / dxdydxdy' (y 

'[-1,1]* 



5.4 The second spatial derivatives 

Let ^ = (£, 77). The real part of £ 2 (£, v) is a correction to e(£, rj) = £77. Since 
d^d v e(^, 77) = 1, we calculate the correction that S 2 gives to that quantity. 

Lemma 5.2 For small q 7^ 

Q2 

lim Re— — S 2 (g ,0) = (2 + 4 log 2) (log |g |) 2 + 0(| log |g | |) (23) 

Proof: By symmetry it suffices to consider q > 0. In general, the second order 
spatial derivatives are 



-d i d j Z 2 ( qo ,q) = z[ i ' j) + Z^ + zi i ' j) 



where 



igo + £ / (i?o + e) 



2 



4 J) = (2(9 l e)(9,e)(-^)(E 1 ) 
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Here e = e(k 2 — k 3 + q) = E 1 and e — E 2 — E 3 — E ± . Denote 

C(«.) = -^E 2 («b,0,0) 

In the xy case, the two integration momenta are denoted by (x, y) and (#', y'), and 
£ 2 = xy, E 3 = x'y', and i?i = e(k 2 - k 3 + q) = + x - x')(r] + y - y'). Since 
% \^= v =o= (y - y')> <V \n= v =o= (x - x 1 ), and d^e = 1, 



C = Ci + C 2 + C 



3 



with 



f \ A*V (\ I 9 (^)fa-i/) A (/l+&23)(/2-/3) 

SI - l^,, \ ^ Z iq +s(x,y,x',y')J (iq +s(x,y,x> ,y>W 



C2 



)) 2 

/_ ii]4 d4 * (-m^o-^)) i^ig&F) 

/_ 11]4 d4 * 2 (-^)<W) Oiw 



Here X = (x, y, a;', y') and d 4 X = dx dy dx' dy', and 

e(x,y,x',y') = xy - x'y' - E x = xy - x'y' - (x - x')(y - y'). 
Using the decomposition given in Appendix B, 

Ci — > - 2 / dxdydx'dy' V f 1 + 2 , Fj ^ ^— - 11 (p,) 

(The limit can be taken under the integral. Decomposing according to the signs 
of x, x' . . ., only j E {1, 2, 3, 4, 9, 10, 11, 12} can contribute. Using the symmetry 
(31), one obtains the above. The minus sign arises from the combination of Fermi 
functions fp, as discussed around (16).) We write Ci = — (Cn + £12)- The term 
involving the 1 is 

r 4 1 

Cn = 2 / dxdydx'dy' £ - 11 ( Pi ) 

By (32) and (33), 

Cn = 8 / dxdydx'dy' Re - \ 1 ( Pl ) . 
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This is the same (up to a sign) as (17), (18), hence £ii = 41og2 (log |<7o I ) 2 + 
O(llogkll). 

The term £ 12 involving the Fj gives another contribution, 

C12 = 2 / dxdydx'dy' £ 2 g 11 

The summand for j = 2 gives the same integral as that for j = 1 because the 
integrand is related by the exchange x <-> y and x' <-> y'. Ditto for j = 4 and 
j = 3. Thus 

C12 = 8 / dxdydx'dy' £ g- 11 
i[o,i] 4 jti^a ^° + e i> 

Because F x = —F 3 and e x = —e 3 , and because p 1 p 3 , 



C 12 = 16Re / dy / dy' / dx' / dx (x X ' )(y + y/) 



Thus 

-1 r l /•! rx' 

jo jo Jo ~" (ig + x'(y + y')+y'(x'-x)) 3 

Let 6 = x'(y + y'), change integration variables from x to u = x' — x e [0, x'] 
and use 

r' dM « = ^£ 

7 (ig + & + y'M) 3 2(ig + fe + xy) 2 (ig + 6) 
Renaming to z = x', we have 

1 " fl (y + y')^ 



C " = - 8 l iy l iy 'l dzRe (^(y + y<))(i* + ;(y + 2y<))* 
The bound 

(y + y') z 2 



< 



q 2 + z 2 {y + 2 yy 



(ig + ^(y + y')) (ig + <Y + 2y')) 2 
for the integrand implies that 

"1 /•! 1 / /■„ , o Tr /\2 



(y + V) 



% 2 
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which shows that |£ 12 | < const [log |g |] 2 for small \q \. 
We now calculate the coefficient of the log 2 . First rewrite 

(j + j')z 2 1 z 2 



(iq + z(y + y')) (iq + z(y + 2y')) 2 (y + 2y') 2 (z + \A){z + LB) 2 

with 

A = ^L- B= qo 



y + y' "" y + 2y' 

Partial fractions give 

£ 2 « /3 . b 2 

+ 1 



(,2 + L4)(,2 + LB) 2 z + iA z + iB (z + LB) 2 
with 

A 2 /y + 2y'^ 2 



(-B - A) 2 V y' 

£(B-2A) , /y + 2y' 



P = ~Tn 7vr = l 



(s-A) 2 V y' 

B 2 „y + y' 



2 



5 



A — b y 

a, (3 and 6 2 are real, so we need 

" x dz 1 



and 

Re io 2 



o (,2 + LB) 2 5(1 + B 2 ) 



With this, we have 



1 pi 



C12 = -4 / dy / dy' 1 



o (y + 2y') 2 

aln(l + A- 2 )+/31n(l + fi- 2 ) + 



26o 



5(1 + fi 2 ) 
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Collecting terms and renaming y' — rj gives 

C12 = -4 / dy f d V 
Jo Jo 



T? 2 ? 2 + (y + 277) 2 



+2 



+ 



y + ?? i 

77 q 2 + {j + 2 V y 



(y + 2 V y 



In 1 + 



(y + 277) ; 



Although the first summands individually contain nonintegrable singularities at 
rj — 0, these singularities cancel in the sum. A convenient way to implement this 
is to use that 



1 ln Qo + (y + v) 2 

T? 2 g2 + (y + 2r?) 2 



+ 2- 



V Qo + (j + 2 v) 2 
d_ 

drj 



_i ln % 2 + (y + ^) 2 



+- 



r/ 



?o + (y + 2r ?) 2 . 

2(y + 77) 2(y + 4r/) 



<? 2 + (y + v) 2 Qo + (y + 2t7) s 



(24) 



Moreover 



1 



(y + 2r/) 2 



In 1 + 



<i 



1 _d_ 

2 ^77 



y + 2?7 



In 1 + 



Thus 

Cl2 = 



dy 



.1 ln % 2 + (y + ^) 2 



(y + 27y) ; 



1 1 



+ 



g 2 + (y + 2r / ) 2 



77 g 2 + (y + 2 V ) 2 2y + 27] 



In 1 + 



(y + 2r?) ; 
% 2 



J7=l 
??=0 



—4 / dy / d?7 



n 



y + v 



y + 2 V 



ql + (y + v) 2 il + (y + 2^ 



Evaluation at 77 = 1 gives one bounded term and one term of order log \ q \. The 
terms at 77 = give 



./o 



dy 



2y 



.2 ^ , r2 2y 



9d + y 



- — In I 1 



% 2 
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The first summand integrates to 0(ln |g |). By (21), the second term gives 

2 (In | g 1) 2 + bounded 

The remaining integrals are 

dy / dr] 



8 / dy / d V 



ql + ij + v? 

since the integrand is bounded by 1, and 



% 2 + (y + v) 2 <7o + (y + 2??) 2 
l 

< 8 



/dy f 

Jo Jo 



drj - 
V 



y 



y 



ql + (y + v) 2 ql + (y + ?v) 2 



The integrand is bounded by a constant times g a +y 2 +T; 2 so the integral is of order 

log|g |. Thus C12 = 2(log|g |) 2 + bounded terms, so that = -((u + ( 12 ) = 
—(4 log 2 + 2) (log |go|) 2 + less singular terms. 

We now show that ( 2 and £3 vanish in the limit (5 — > 00. Consider ( 2 first. Let 
G P (E!) = -£k%(£i) - ^(£1). (Note that G P (E) = pG^pE).) At £ = rj = 0, 
i?i = (rr — — y') is invariant under the exchange (x, y) <-> (#', y'), so 



C: 



d 4 X Gp{E 1 )f p {xy) ( : 
i,i] 4 



1 



1 



Vig + e(a:, y, a:', y') ig + e(a:', y', x, y) 
Because e(x, y, x', y') = xy' + x'y — 2x'y', 

e(x', y', x, y) - e(x, y, x', y') = 2(x'y' - xy) = 2[x(y' - y) + (x' - x)y'\ 
Hence 

x(y' -y) + (x' - x)y' 



(2 

-1,1] 4 
Consider first 



d 4 X Gp(E 1 ) fp(xy) 



Ti = 2 



-i,i]' 



d 4 X G^) ^(zy) - 



(ig + e(ar, y, a;', y')) (ig + e(z', y', x, y)) 
x(y' - y) 



(ig + e(ar, y, a;', y')) (ig + e(a/, y', x, y)) 
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Because [— 1 , l] 4 and the integrand are invariant under the reflection R which maps 

(x, y, x', y') to (-x, -y, -x', -y') and 

R{(x,y,x',y') \ y > y'} = {(x,y,x',y') | y < y'}, 
we can put in a factor 2 1 (y > y') . Changing variables from x' to E\, so that 

Ei 

x = x 



y-y' 



gives 



Ti = 4 / dxdydy' 11 (y - y' > 0) (-x) f p {xy) 



(x-i)( y -y>) 1 (m + e(x, y, x', y')) (iq + s(x', y', x, y)) 
Change variables one last time, from E 1 to u = @Ei, to get 

Ti = 4 / dxdydy' / du 
J[-i,i] 3 Jm. 

l(y-y'>0) 11 (P(x -l){y-y')<u< (5{x + l)(y- y')) 

Gi(u) 



(~x)fp(xy) 



(iq + e(x, y, x', y')) (iq + e(x', y', x, y)) 



where /3 — > oo. The integrand is bounded in magnitude 

by the L 1 function 

(x,y,y',u) i-> Gi(u)-2 
and it converges almost everywhere (namely, for x ^ ±1, xy ^ 0) to 

1 (y - y' > 0) x 9N) ^^ ( ^, )2 K-#;(,) - 

By dominated convergence, the limit (3 — > oo can be taken under the integral. The 
limiting range of integration for u is R. By the fundamental theorem of calculus 

/ duG 1 (u)= [ du^[-u5 1 (u)] = 0, 
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so Ti vanishes as (3 — > oo. The calculation of the limit (3 — > oo of 

7 1 — 2 / d 4 X G (E ) f (xv) — — 

is similar and gives as well. Thus ( 2 — > as (3 — > oo. By the same arguments 

where £ = y, x', y') and e = y' , x, y). After the same limiting argument 
as before, the w-integral is now 

/ u8\{u) du = 
Jr 

because the integrand is odd. Thus ( 3 — » as (3 — > 00, too. ■ 



Lemma 5.3 77ze rea/ part o/7/ze second derivative of the self-energy with respect 
to £ grows a? most logarithmically as go - ► 0: ?/zere are constants A and B such 
that for all |go| < 1 



< A + B lot 



9o 



77ze same /zo/ds ov symmetry for the second derivative with respect to rj. 
Proof: We need to bound 



9, 



2 / (MEh) + b p {E 2 - E 3 )) (f p (E 2 ) - f p (E 3 )) 



igo + ^ 



at £ = rj = where all ^-dependence is in E x — (£ + x — x')(?7 + y — y') and in 
e = E 2 — E 3 — Ei = xy' + x'y — 2x'y'. We proceed as for the mixed derivative 
d 2 /d£dr/, but now some terms are different because = 0. We obtain 



d 2 S 2 
d£ 2 



(g , 0, 0) = Xx + X 2 + X 3 
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with 



Xi = 2((y-y 
X 2 = 



A2 + ^ - e 3 )) {f p {E 2 ) - ^(^3)) 



(i<?o + ef 



l\2 



= 2 



(y-y') 
iqo + £ 

MEi) (y-y') 2 



We first calculate the zero-temperature limit of X 1 . Using the notations of Ap- 
pendix B, 

XS 0) = lim Xx = -4 / dxdx'dydy' V D " 11 (p n ) 

We use that .D 2 = £> 2 and e 3 = —e±, to combine n = 1 and n = 3 in one term, 
and n = 2 and n = 4 in another, so that xf } = xf J + X^ with 



X 
X 



S = -8i f dy f dy' (y' + y) 2 f dx' f dx Im (i9()+x , y+y , (2x ,_ x 
S = ~8i ^ [ dy (y' - y) 2 £ dx' £ dx Im (igo+xy , +x , (2y , 



)) 3 
y)F 



Thus Xi does not contribute to Re <9|£ 2 . 

We now show that X 3 is imaginary in the limit (5 — > 00, because after a rewrit- 
ing of terms, the difference f/3(E 2 ) — fi3(E 3 ) effectively implies taking the imag- 
inary part. Because (y — y') 2 Sp(E 1 ) is invariant under (x, y) <-> (x',y') (and 
denoting e(x, y, x', y') = e(x', y', x, y)) 



h = -ij^Xiy-y'fSpmfpixy) 



1 



(iq + e) 2 (iq + e) 2 
= -A j dy j dy' l{y-y' >U) J dx f p {xy) (y - y') 

[ dx' (y - y')^(^i) L 1 V 2 - c 
J l{iq + e) 2 {iq + e) 2 _ 

For the second equality, we used invariance under the reflection (x, y, x', y') — > 
(— x, —y, —x', —y'). The convergence argument used in the analysis of the 7\ 
contribution to d| £ 2 can be summarized in the following Lemma. 
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Lemma 5.4 Let [3 > and F : \J3 , oo) x [—1, l] 4 —> Cbe bounded, and 
lim j^r y y') = f(x,y,y') 

a.e. in (x, y, u, y'). Let E 1 — {x — x')(y — y'). Then 



Jim f F((3,X)5 p (E 1 )(y-y') l{y > y') d 4 X 
dy J dy' l(y> y') J dx f(x, y, y') 



Applying Lemma 5.4, and using that e{x, y, x, y') = x(y — y') = —e{x, y, x'y'), 
we get 

i y i 
jirnX 3 = -A J dy J dy' (y - y') J dx 0i (-xy) l(x, y, y') 



-l -l 



-l 



where 0i(i) = lim fp(—x) is the Heaviside function, except that 0i (0) = 1/2, 

2 /3^oo ' 2 



z(x,y,y') = 



l 



l 



= 2ilm 



(iq + x(y - y')) 2 (iq - x(y - y')) 2 
1 



(iq + x(y -y')) 2 ' 

Thus, in the limit (3 — > oo, X 3 does not contribute to the real part of <9|£ 2 either. 

It remains to bound X 2 . Here we integrate by parts, to remove the derivative from 
the approximate delta function, and then take (3 — > oo. 
We first do the standard rewriting 



X 2 = 2 J d 4 Xt(y-y'>0) {y-y'f 
fp{xy)5'p{{x-x'){y-y')) 



1 



iq + e iq + e 



We integrate by parts in x, using that (y — — x')(y — y')) = ^^{{x 

x')(y — y 1 )). This gives the boundary term 

B w = 2^x Jdyjdy't(y-y'>0) J dx' f p {xy) 

1 1 



(y - y')8p((x - x')(y - y')) 



iqo + e iq + e 
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and two integral terms, namely 

h = 2 j dy j dy'l{y-y'>Q) J dx J dx' yS^xy) 
(y- y >)5p{{x-x'){y-y')) 

and 



ig + £ i<?o + £ 



I 2 = 2jdyjdy'l(y-y'>0) J dx J dx' f p (xy) 



(y- y >)5p((x-x')(y-y')) 



y' y' - 2y 

_{iq + e) 2 {iq + ef 



An obvious variant of Lemma 5.4, where x is summed over ±1 instead of inte- 
grated, applies to the boundary term and gives 



i y 



B = lim B 1/f) = 2 V x I dy I dy' Qi(-xy) 



-l -l 
1 



1 



_iq + x(y - y') iq - x(y - y') 
Thus B is real. Because 6i (y) + 0i (— y) = 1, we get 

i y 



B = ijdyj dy' 



y-y 



<i + {y- y') 



With z = y — y', we have 

1 1+2/ 

B = ijdyjdz 

-1 o 

i 



ql + z* 



= 4 J dy i ln(g 2 



-1 

= O(log|g |). 



l+y 


2 / drj In ( 1 + ^ 
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In Ii, we change variables from (x' : x) to (u, v) where u = (3(x — x')(y — y') and 
v = (3x\y\ and get 



P\v\ 



h = 2 J dy J dy' sgn(v) J dv S^v) J du 6i(u) 



-1 -1 



-P\v\ 



1 



1 



1^0 + £ (j\y~\ ) Vi ~fi\y] ~ l3(y U -y>) ' y') ^° + 



P\v\ /%-2/') ' y ' 



The integral converges in the limit (3 — > oo by dominated convergence. The last 
factor vanishes in that limit, so 1\ — > as /3 — > oo. 

Finally, Lemma 5.4 implies that = lim^oo I 2 exists and equals 

i y i 



if = 2jdyjdy' J dxQii-xy) 



-l -l 



J/'-2y 



(i^o + ^w-w')) 2 (m-x(y-y')) 2 
In the real part, the terms with y' in the numerator cancel, so that 

i y i 



Re lf ] = 2 J dy J dy' J dx Gi(-xy) Re — 



2y 



-i -i 
i i 



= 2 dy dx Qi(-xy) Re 



1+2/ 



dz- 



x(y - y')) 2 
2y 



i i 



= —4 J y dy J dxQi(-xy) 



= -4 



-l -i 

1 1+2/ 

jyiy J 

"I -(1+2/) 
1 1+2/ 

= -4 y y dy 

1-2/ 



dv 6! 



o " (ig - xz) 2 

5 2 +a . 2(1 + y)2 

r/y \ 1 



2 \ 1+!// ?o+" 2 



dv 



ql + v 2 
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The contribution from the region v > 1 is obviously bounded. The remaining 
integral is 



11 11 
dv f dv 



fyiyj 



ql + v 2 J Qo+ v 2 



y dy 



l-j/ 1-v 

1 



5.5 One-loop integrals for the xy case 

For the discussion in Section 6, it is useful to calculate the lowest order contribu- 
tions to the four-point function, the so-called bubble integrals. Again we restrict 
to the xy type singularity. Because the fermionic bubble integrals are not contin- 
uous at zero temperature, it is best to calculate them by setting q = first, then 
letting the spatial part q tend to 0, all at a fixed inverse temperature f3, and then 
calculate the asymptotics as (3 — > oo. 

s + q — «-^- < pp*- t + q q — s — »^ > y^- q — t 

s — > £ s — ^ > £ 

particle-hole bubble particle-particle bubble 

5.5.1 The particle-hole bubble 

Write x' = x + £, y' = y + rj. The bubble is 



dxdy - V : ■ ; 

, -1,113 p luj - x y Aqo + 



) - x'y 1 



L 



dxdy f ^ XV) ~ f ^ X ' Vl) 



[-i,i]2 igo + xy- x'y 



= [ dt [ dxdy . Xy Xy - , {-8p){txy + (1 - t)x'y') 
Jo J[-i,i] r ' : '" ~ ' " 



12 iq + xy - x'y' 



(25) 
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Obviously, at q ^ 0, B ph (q , £, rj) — > as (£, 77) — > 0. So set g = 0, i.e. consider 

<h= I™ S ph (0,e,77). 

Lemma 5.5 77?e /arge /3 asymptotics of is 

B^ = -2\n(3 + 2K + 0(e- f3 ) 



where K = j^rulnu. 
Proof: By (25), 



= / dxdy {-5p)(xy) = 4 / dxdy (-^(xy) 
J[-i,il 2 ^fo,il 2 



,i] 2 ^[0,1] 
f 1 da; j px du | ^ dw /3 

io a? io 4 cosh 2 £ Jo 4 cosh 2 f it 



2 </ u wu " 2 

'' x d// u 



-2\n(3 I" dU 2 +2K- [ 
H Jo 2 cosh 2 2 L 



hi—. 

cosh 2 s /3 



2 J/3 cosn 2 

The last integral is exponentially small in f3 because of the decay of 1/ cosh 2 . 

5.5.2 The particle-particle bubble 

This time write x' = x — £, y' = y — rj. The bubble is 

/I y y 
dxdy - V : r — 
t -i,i] 2 p „ - 1(50 - - ^ y 

fp(-xy) - fp(x'y') 



dxdy 



i,i]2 " -i<?o + xy + x'y' 

1 , (P 

dxdy tanh — xy 

(ta)-o Jr 1 112 -ig + 2xy V 2 



Again, we set q = and keep /3 < 00. Then 

Bg? = 5 PP (0,0,0) 



pp 





dxdy tanh ( —xy) = 2 [ dxdy — tanh 

i,i] 2 2x V V 2 / i[o,i] 2 xy 



2 ^ 
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Jo x Jo 



— tanh-E 
E 2 

1 



- 2 ; ' ,!/•; — tanh -/•; 

£ 2 



-(ln£) 2 tanh^E 



o 



+ / dEQnE) 



P 



2 cosh 2 §E 



1 



/3 / cosh 2 v 

2 



Jo \ (3 ) cosh 2 -y 



2.. I dv ( ln 5/ _u2 



/3 / cosh v 



Thus we have 



Lemma 5.6 



= (ln/5) 2 - 2inn/3 + AT' + O (e~ p ) 
where K = ( n °° ^Ldv and K> = f n °° ^^dv. 

JD cosh ti -"U cosh d 



6 Interpretation 

Let us discuss these results a bit more informally. The above calculations for 
the xy case can be summarized as follows. Evidently, the derivatives we were 
looking at diverge in the limit q — > at zero temperature. To proceed, we discuss 
positive temperatures and replace q by tt//3, and thus translate everything into 
/3-dependent quantities. We have proven that to all orders r in A, 

|VE r | < const (26) 

(where the constant depends on the order r in A). In the model computations of 
the last section, we have seen that to second order in the coupling constant A 

-4 In 2 (Aln/3) 2 
(2 + 4 In 2) (Aln/3) 2 

0(A 2 La/3) 
-2 Am/3 

A(ln/3) 2 -|(Aln/3) 2 



Im<9 X 
Re 



a 2 * 

d£dri' 
a2 



Re 

5 P h(0) 

^pp(O) 
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We have redefined B ph and B pp to include the appropriate coupling constant de- 
pendence. The line for B pp includes second order corrections to the supercon- 
ducting vertex, where the coefficient given by the loop integral is a. If A is nega- 
tive (attractive interaction), the superconducting instability driven by the A(ln (3) 2 
term is always strongest, but if A > (repulsive bare interaction), the \(\n(3) 2 
term suppresses the leading order Cooper pair interaction, and the higher order 
term proportional to a is only of order (A In (5) 2 . In this case, all terms that can 
drive instabilities are linear or quadratic in A In (5. A first attempt to weigh the 
relative strength of these divergences is to look at the prefactors. Here it seems 
that the second derivative gets the largest contribution. The asymmetry between 
this logarithmic divergence and the boundedness of the gradient is striking. 

In the following we discuss the physical significance of the above findings. 
This discussion is not rigorous, but it reveals some interesting possibilities that 
should be studied further. 

We first note that the above results were achieved in renormalized perturbation 
theory, that is, the above all-order results are in the context of an expansion where 
a counterterm is used to fix the Fermi surface. (The second order explicit calcu- 
lations assume only that the first order corrections can be taken into account by a 
shift in /j, which is true because the interaction is local.) Using a counterterm is 
the only way to get an all-order expansion that is well-defined in the limit of zero 
temperature. A scheme where the scale decomposition adapts to the Fermi surface 
movement was outlined in [15] and developed mathematically in [16] and [17]; 
but in any such scheme the expansions have to be done iteratively and cannot be 
cast in the form of a single renormalized expansion, because the singularity moves 
in every adjustment of the Fermi surface [15]. The meaning of the counterterm 
was explained in detail in [10, 11, 12, 13]. In brief, the dispersion relation we 
are using in our propagators is not the bare one, but the renormalized one, whose 
zero level set is the Fermi surface of the interacting system. Thus, we have in fact 
made the assumption that 

the Fermi surface of the interacting system has the properties H1-H6. 
In particular it contains singular points and these singular points are 
nondegenerate. 

We shall first discuss our results under this assumption and then speculate about 
how commonly it will be valid. In the case of a nonsingular, strictly convex, 
curved Fermi surface, there was a similar assumption, which was, however, math- 
ematically justified by our proof of an inversion theorem [ 1 3] that gives a bijective 
relation between the free and interacting dispersion relation. 
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6.1 Asymmetry and Fermi velocity suppression 

First note that the regularized (discrete-time) functional integral for many-body 
systems has a symmetry that allows one to make an arbitrary nonzero rescaling of 
the field variables. This is based on the behaviour of the measure under ipi — > g i ^ i 
and — > The result is 

J fj djityi e- A{ ^ = (j[ gfg)j ]J d^d^ e"^" 1 ^- V) (27 ) 

Source terms get a similar rescaling. This can, of course, be used to remove 
a factor Z^ 1 from the quadratic term in the fields, but the factor Z (see (13)) 
will then reappear in the interaction and source terms. Note that Z depends on 
momentum. In the limit of very small Z, some terms may get greatly enhanced or 
suppressed. 

The Fermi velocity is defined as v F (p) = ZVY* on the Fermi surface. If one 
extrapolates the above formula for <9 S2 to (3 — > oo, the Z factor, defined in (13), 
becomes zero at the Van Hove points. 

There is a crucial difference between the one- and two-dimensional cases. In 
dimension one, both <9 22 and VS2 behave like A 2 log (3, so that, after extract- 
ing the field strength renormalization, the Fermi velocity retains its original value. 
In our two-dimensional situation, however, VS remains bounded, and thus the 
Fermi velocity gets suppressed in a neighbourhood of the Van Hove singularity 
because it contains a factor of Z. Because the time derivative term in the action is 
Z^ipikoip and k is an odd multiple of the temperature T = (3~ l , one can also in- 
terpret Z{k)~ l T = T(k) heuristically as a "momentum-dependent temperature" 
that varies over the Fermi surface and that increases as one approaches the Van 
Hove points ("hot spots"). This behaviour is illustrated in Figure 1. 

6.2 Inversion problem 

There is, however, also a more basic problem that is exhibited by these results. 
It arises when one starts questioning the assumption that the interacting Fermi 
surface contains singularities. The second derivative of the self-energy in spatial 
directions is divergent at zero temperature. Thus even the lowest nontrivial cor- 
rection may change the structure of e significantly. This is related to the inversion 
problem, which was solved for strictly convex Fermi surfaces in [13]. 
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-3 

-3 -a -i 01 i o 

Figure 1 : Comparison of the derivative of the Fermi function, where on the right, 
the inverse temperature (5 is scaled by an angle-dependent factor that vanishes 
near the points (ir, 0) and (0, tt). The dark regions correspond to small derivatives, 
the bright ones to large derivatives. 

In all of the following discussion, we concentrate on two dimensions and as- 
sume that by adjusting fi, we can arrange things such that the interacting Fermi 
surface still contains a point where the gradient of the dispersion relation vanishes. 

In the general theory of expansions for many-fermion systems, the most sin- 
gular terms are those created by the movement of the Fermi surface. Counterterms 
need to be used to make perturbation theory well-defined, or an adaptive scale 
decomposition has to be chosen. To show that the model with counterterms cor- 
responds to a bare model in some desired class, one needs to prove an inversion 
theorem, which requires regularity estimates. To get an idea, it is instructive to 
consider a neighbourhood of a regular point on the Fermi surface. By an affine 
transformation to coordinates (u, v) in momentum space, i.e. shifting the origin 
to that point and rotating so that the tangent plane to the Fermi surface is given by 
u = 0, the function e can be transformed to a function e(u, v ) = u + |f 2 + . . ., 
where k denotes the curvature and . . . the higher order terms. (By an additional 
change of coordinates, the error terms can be made to vanish, but this is not im- 
portant here.) Let So(p) be the self-energy at frequency zero (or ±7r//3). Then 

(e + £ )~(u, v) = u + -v 2 + Xa(u, v). 

In order for the correction Act not to overwhelm the zeroth order contribution e, 
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we need d u a to be bounded and d 2 a to be bounded. Then the correction is small 
when A is small. Note that different regularity is needed in normal and tangential 
directions. Indeed, even in the absence of Van Hove singularities, the function a 
is not twice differentiable in u at T = 0, but it is C 2 in v if e has the same property 
[11]. 

Thus, the mere divergence of a derivative does not tell us about a potential 
problem with renormalization; the relevant question is whether the correction is 
larger than the zeroth order term. In the strictly convex, curved case, there is 
no problem. However, in two dimensions, the normal form for the Van Hove 
singularity, e(uv) = uv gets changed significantly because 

(2 = £r v Hu,v) + a(u,v)}\ u=v=0 (28) 

diverges at zero temperature. To leading order in q , u and v , our result for the 
inverse full propagator to second order takes the form £2 iq — ( 2 uv with 

C 2 = l + tf 2 (Aln/3) 2 
C 2 = l + ^Am/T) 2 . 

Because ( 2 diverges in the zero-temperature limit, it is not even in second order 
consistent to assume that the type of (u, v) = (0, 0) as a critical point of e and 
E = e + cr is the same. One can attempt to fix this problem by rescaling the 
fields by Q 1 / 2 . Then the iq term gets rescaled similarly because ( depends on 
the same combination of A and In (3 as (. It should, however, be noted that, as in 
the above-discussed suppression of the Fermi velocity, a corresponding rescaling 
of all interaction and source terms also occurs^ In particular, a, which gives the 
Fermi surface shift, itself gets replaced by <rC _1 . Because a remains bounded 
(only its derivatives diverge), this would imply that the Fermi surface shift gets 
scaled down in the rescaling transformation, indicating a "pinning" of the Fermi 
surface at the Van Hove points. 

In our case, $ 2 = 4 In 2 and $ 2 = 2 + 4 In 2. With these numerical values, 
(2 is larger than £ 2 , so ^ 2 appears as the natural factor for the rescaling. In the 
Hubbard model case, the transformations leading to the normal form depend on 
the parameter 6 in (1), so that one can expect ( and £ to depend on 6. The study 
of these dependencies, as well as the interplay between the two critical points that 
contribute, is left to future work. 

One point of criticism of the Van Hove scenario has always been that it ap- 
pears nongeneric because the logarithmic singularity gets weakened very fast as 
one moves away from the Van Hove densities. If the above speculations can be 
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substantiated by careful studies, the singular Fermi surface scenario may turn out 
to be much more natural than one would naively assume. Moreover, there is a nat- 
ural way how extended Van Hove singularities may arise by interaction effects. 



A Interval Lemma 

The following standard result is included for the convenience of the reader. 

Lemma A.l Let e, r] be strictly positive real numbers and k be a strictly positive 
integer. Let I C K be an interval (not necessarily compact) and f a C k function 
on I obeying 

\f {k \x)\ > 77 for all xe I 

Then 

Vol{xeI\\f(x)\<e}<2 k+ \^ k . 
Proof: Denote a = (fj) 1 ^. In terms a, we must show 

\f {h) (x)\ > ^ for all x e / =>■ Vol{ x e I \ \f{x)\ < e } < 2 k+1 a 

Define c k inductively by c 1 = 2 and c k = 2 + 2c k -i. Because b k = 2~ fc c fc obeys 
bi — 1 and b k = 2~ k+1 + b k _ 1 we have b k < 2 and hence c k < 2 k+1 . We shall 
prove 

|/ {fc) (^)l>^ for all x El =}► Vo\{xeI\\f(x)\<e}< Ck a 
by induction on k. 

Suppose that k — 1 and let x and y be any two points in{xe/||/(x)|<e}. 
Then 

\x-v\ = ^ I fix) - f(v)\ = < 2e 

1 yi \f(x)-f(y)\ u {X) ny)l l/'(C)l - l/'(C)l 

for some C e /. As |/'(C)I > i we have \x-y\ < 2a. Thus { x e 1 1 \ f(x)\ < e } 
is contained in an interval of length at most 2a as desired. 

Now suppose that the induction hypothesis is satisfied for k — 1 and that 
|/( fe )(:r)| > ^ on I. As in the last paragraph the set { x G / | \f^ k ~ l \x)\ < } 
is contained in a subinterval J of / of length at most 2a. Then I is the union of 
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Iq and at most two other intervals i+, /_ on which \f^ k ^(x)\ > ^f=t. By the 
inductive hypothesis 

Vol{xe/| \f(x)\ <e} < Vol(/ ) + X! Vol {^e^l 1/0*01 <e} 

< 2a + 2cfc_iQ; = c^a 



B Signs etc. 

The restrictions (16) are summarized in the following table 

E 2 = xy E 3 = x'y' E 1 = (x - x') (y - y 1 ) 



+ 

+ + 



(-1) (29) 
(-1) 



In both cases, the product of indicator functions resulting from the limit of Fermi 
functions is —1. 

Let R be the reflection at zero, 

R(x, y, x', y') = {-x, -y, -x', -y'). (30) 

The function 

£ = e(x, y, x', y') = xy' + x'y - 2x'y' 

satisfies 

s(x, y, x', y') = e(-x, -y, -x', -y'). 
The function D{x, y, x', y') — y — y' satisfies 

D(R(x,y,x',y')) = -D(x,y,x',y') 

The function F( 

')(y — y') is invariant under R. 
In the following we list all cases of signs for x, x', y and y', together with 
e,D,F written as functions of 

x = \x\,y= \y\,x = \x'\,y' = \y'\. 

to be able to restrict the integrals to [0, 1] whenever this is convenient (by trans- 
forming to x, . . . , y' as integration variables), and to exhibit some important sign 
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changes. In the last column, we list the condition p n obtained from the restriction 



on the sign of ( 

n xy x' y' 


x-x')(y-y') in (29). 
£n — xy' + x'y — Ix'y' 


D n = y-y' 




Pn 


1 + + +- 

2 + + -+ 

3 H h+ 

4 + 


ei = x'y+(2x'-x)y' 
e 2 = xy' + (2y' - y)x' 
£ 3 = xy' - (2y' + y)x' 
£ 4 = x'y - (2x' + x)y' 


Di = y + y' 
D 2 = j-y' 
D 3 = -(y + y') 
D, = -(j-j') 


(x-x')(y + y') 
(x + x')(y-y') 
-(x-x')(y + y') 
-(x + x')(y-y') 


x < x' 

y < y' 

/ 

X < X 

y < y' 


5 + + ++ 

6 + + 

7 H 1 — 

8 +--+ 


£ 5 = xy' - x'(2y' - y) 
£ 6 = -(xy'+x'(2y' + y)) 
e 7 = -(xy' -x'(2y'-y)) 
e 8 = xy' + x'(2y' + y) 








9 + 

10 --+- 

11 - + 

12 - + ++ 


£i 

£2 
£3 
£4 


-D\ 

-D 2 
-D 3 
-D 4 




pi 

P2 
P3 
Pi 


13 

14 --++ 

15 - + -+ 

16 - + +- 


£5 
ee 
£7 

£8 









Cases 1-4 and 9-12 obey the restrictions (29). Cases 5-8 and 13-16 do not because 
there, the signs of xy and x'y' are the same. They are used to discuss some terms 
at finite (3. Since the first two restrictions are not satisfied, the column for the last 
restriction, p, is left empty in these cases. Case n + 8 is obtained from n by the 
reflection R. Thus 



£n+8 — £n, Ai+8 — —D n , F n+8 — F n , p n+8 — p n (31) 

Moreover, 

Si = —e 3 and £ 2 = —£4 (32) 

and 

£ 2 (y,x,y',x') = £i(x,y,x',y') and p 2 (y, x, y', x') & Pl (x, y, x', y'). (33) 
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